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APPENDIX  i 

TRAJECTORY  ANALYSIS  AND  SYNTHESIS 


Introduction 


The  mission  of  the  scoring  methods  study  is  to  investigate  those 
phenomena  which  may  have  potential  in  the  development  of  new  scoring 
techniques,  and,  ultimately,  to  develop  concepts  for  new  methods  of 
determining  (scoring)  the  relative  trajectory  of  a,  munition  with  respect 
to  its  target.  One  of  the  problems  associated  with  this  study  is  that  of 
determining  how  (or  whether)  useful  trajectory  information  can  be 
generated  through  full  exploitation  of  existing  mathematical  techniques. 
More  specifically,  what  are  the  characteristics  of  the  relative  tra- 
j-actories  and  how  can  they  be  used  advantageously  in  the  development 
of  a  trajectory  scorer?  What  are  the  limitations  of  the  assumptions 
that  it  would  be  convenient  to  make? 

This  paper  deals  with  the  use  of  position,  velocity,  and  accelera¬ 
tion  data  to  give  trajectory  information  or  improvements  over  meas¬ 
ured  trajectories.  The  constant-relative-velocity  assumption  is 
investigated  and  the  velocity  and  acceleration  characteristics  of  a 
rocket  are  developed  with  atmospheric  drag  and  gravity  effects  neg¬ 
lected.  The  problem  of  non-zero  relative  acceleration  is  considered; 
and  a  particularly  useful  method  for  smoothing,  interpolating,  extrap¬ 
olating  and  (if  desired)  integrating  or  differentiating  trajectory  data 
is  discussed  briefly.  The  topics  just  mentioned  are  taken  up  in  the 
order  named  after  a  preliminary  discussion  of  trajectory  synthesis. 

Synthesis  of  Trajectories 


The  term  "trajectory  synthesis"  will  be  used  to  mean  the 
techniques  and  procedures  by  which  observed  data  are  used  to  recon¬ 
struct  the  path  of  an  object  over  some  prescribed  interval  of  observa¬ 
tion.  This  definition  includes  any  operations  needed  to  obtain  values 
of  the  trajectory  variables  from  the  reduced  data  and  any  data-improve- 
ment  activities  such  as  smoothing,  curve  fitting,  interpolation,  and 
extrapolation. 

For  scoring  applications,  the  desired  result  is  the  relative  tra¬ 
jectory  of  the  weapon  with  respect  to  its  target.  This  can  be  obtained 
directly  by  observing  the  motion  of  the  missile  from  a  station  or 
stations  on  the  target,  and  vice  versa,  or  by  making  some  observations 
from  each  vehicle.  The  indirect  procedure  is  to  measure  the  tra¬ 
jectories  of  the  target  and  the  weapon  with  respect  to  an  external  ref¬ 
erence  such  as  a  coordinate  system  attached  to  the  earth.  The  data 
needed  to  establish  such  trajectories  can  be  obtained  from  one  or  more 
stations  external  to  the  observed  objects  or  by  use  of  inertial 
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instruments  (accelerometers  or  velocimeters)  attached  to  the  objects; 
in  the  latter  case,  additional  obsei’vations  will  be  needed  to  establish 
the  initial  conditions.  The  data  taken  may  include  signals  which  are 
directly  related  to  angles,  distances,  or  linear  and  angular  velocities 
and  accelerations  as  functions  of  time.  Frequently  the  data  are  not 
obtained  continuously  but  at  discrete  times  which  are  usually  and 
preferably  equally  spaced;  the  data  taken  may  cover  the  entire  tra¬ 
jectory  of  the  weapon  in  the  vicinity  of  the  target  or  only  portions  of  it. 

A  knowledge  of  the  behavior  of  the  trajectory  variables  including 
their  ranges  and  rates  of  variation  is  essential  to  selecting  a-method  of 
scoring  and  the  instruments  and  procedures  to  be  used.  The  coverage 
required,  the  accuracy  and  response  time  needed  in  the  instruments, 
and  even  the  data  reduction  and  processing  procedures  are  all  functions 
of  these  trajectory  characteristics.  The  complexity  of  the  equipment, 
the  amount  of  effort  required  in  the  processing,  and  of  course  the  cost 
of  the  whole  operation  can  be  considerably  reduced  through  knowledge 
of  the  simplifying  assumptions  which  are  acceptable  in  a  given 
application. 

The  purpose  of  this  paper  is  to  demonstrate  the  probable  charac¬ 
teristics  of  the  several  types  of  relative  trajectories  which  are  of 
interest  in  this  study.  The  validity  of  the  constant-relative-velocity 
assumption  is  considered  for  all  cases,  and  the  characteristics  of  such 
trajectories  are  discussed.  Since  the  scoring  requirements  for  future 
weapon-target  intercepts  do  not  define  the  tactics  to  be  employed,  a 
discussion  of  the  velocity  and  acceleration  characteristics  of  a  rocket 
(with  gravity  and  aerodynamic  forces  neglected)  is  provided  in  support 
of  conclusions  made  earlier  and  for  use  in  checking  these  conclusions 
when  more  definite  information  about  the  problems  becomes  available. 
Some  remarks  are  made  about  trajectories  with  non-zero  relative 
acceleration  and  about  data  processing  operations  (integration,  interpo¬ 
lation,  etc.  )  for  relative  trajectories.  Errors  in  the  data  due  to 
bending,  own-ship  angular  motions,  and  the  measuring  instruments  can 
obscure  but  obviously  do  not  determine  the  true  nature  of  the  trajectory. 
The  effects  of  these  errors  and  the  procedures  needed  to  ininimize 
these  effects  are  the  subjects  of  separate  studies  made  under  this 
contract.  Although  some  remarks  on  differential  errors  are  made  in 
this  paper,  it  will  usually  be  assumed  that  the  data  are  measured  with¬ 
out  error  with  respect  to  a  non-rotating  (sometimes  referred  to  as 
inertial)  coordinate  system. 

A  trajectory  is  established  when  the  position  of  a  moving  point  is 
known  as  a  function  of  time.  In  general,  three  position  coordinates  are 
determined  to  a  certain  accuracy  (directly  or  indirectly)  at  discrete 
points  in  time.  These  time-position  values  obtained  from  the  original 
data  may  be  smoothed,  extrapolated,  and  interpolated.  The  desired 
functional  relation  may  be  obtained  by  means  of  graphs  or  by  curve  fitting. 
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The  desired  trajectory  data  may  be  obtained  from  successive 
observations  of  position  as  a  function  of  time  or,  if  suitable  initial 
conditions  can  be  obtained,  synthesized  from  either  velocity  or 
acceleration  data  or  both. 

It  is  well  known  that  the  velocity  may  be  obtained  by  the  time 
integral  of  the  acceleration: 

p  t  , 

V(t)  =  V(0)  +  \  V  dt 
0 


Position  may  be  obtained  as  the  integral  of  the  velocity  or  the  second 
integral  of  the  acceleration.  For  any  position  coordinate  X, 

X(t)  =  X(0)  +  C  X  dt 
0 


X(0)  +  r  <  X(0)  +  r  ^  X  dt, 

J  0  ^  n  ^ 


p  dt.. 


=  X(0)  +  tx(0)  +  J  Jj^Xdtjldt. 


If  the  velocities  and  accelerations  are  those  of  the  weapon  rela¬ 
tive  to  the  target,  then  the  corresponding  integrals  for  displacement 
represent  displacements  of  the  weapon  relative  to  the  target;  similarly, 
the  relative  velocities  are  integrals  of  relative  accelerations.  If  some 
point  other  than  the  target,  such  as  the  center  of  the  earth  or  a  point 
on  its  surface,  is  chosen  for  reference,  the  integrals  again  give 
relative  velocities  and  displacements  with  respect  to  the  reference 
coordinate  system. 

If  position  coordinates  are  to  be  obtained  as  functions  of  the  time 
by  integrating  the  velocity,  it  is  necessary  to  know  the  initial  values  of 
the  position  coordinates  and  to  obtain  values  for  the  velocity  components 
as  functions  of  time.  In  the  simplest  case,  treated  in  the  next  section, 
the  relative  velocity  is  constant  within  the  accuracy  requirements,  and 
the  results  of  the  integration  can  be  given  in  closed  form.  For  any 
position  coordinate,  X,  the  result  would  be 

X(t)  =  X  +  X  t 
'  '  o  o 

• 

where  X  is  the  initial  condition  and  X  is  the  constant  velocity  com- 
ponent.  An  important  possible  application  for  such  linear  relations  is 
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the  synthesis  of  relative  trajectories  for  head-on  attacks  against 
an  ICBM  or  satellite.  If  portions  of  the  trajectory  do  not  have  a  data 
record,  linear  extrapolation  or  interpolation  can  be  used  to  compute 
these  missing  portions. 

Similarly,  at  least  two  quantities  are  needed  to  obtain  the  velocity 
by  integrating  the  acceleration.  These  are  sufficient  if  the  acceleration 
is  constant,  or  practically  constant.  These  quantities  are  the  initial 
velocity  and  the  average  acceleration.  Likewise,  a  knowledge  of  the 
initial  position,  initial  velocity,  and  the  average  acceleration  is 
sufficient  to  permit  the  computation  of  position  as  a  function  of  the  time 
by  integrating  the  acceleration. 

Just  as  the  position  could  be  determined  by  integrating  the  velocity, 
so  may  the  velocity  be  determined  by  differentiating  the  position  coordi¬ 
nates  as  functions  of  time.  At  least  two  positions  as  functions  of  the 
time  are  required  for  such  a  differentiation  but  more  are  preferable. 

A  complete  record  of  position  vs  time  thus  permits  the  computation  of 
a  complete  record  of  both  velocity  and  acceleration  vs  time.  If  initial 
conditions  of  position  and  velocity  are  known,  a  complete  record  of 
acceleration  is  sufficient  to  enable  the  computation  of  complete  records 
of  velocity  and  position  as  functions  of  time.  In  general  terms,  if  suffi¬ 
cient  initial  conditions  are  known,  then  a  record  of  position,  velocity, 
or  acceleration  permits  the  computation  of  the  remaining  quantities  in 
this  group.  This  procedure  may  be  required  in  computing  the  tra¬ 
jectory  or  portions  of  it  in  scoring  applications. 

A  position  coordinate  obtained  by  integrating  a  velocity  is  subject 
to  an  error  which  is  the  integral  of  a  velocity  error.  If  the  permissible 
error  in  such  a  position  coordinate  is  f  and  the  maximum  time  period 
for  the  integration  is  At,  then  the  maximum  permissible  constant  error 
in  velocity  is  f/At.  Similarly,  if  position  is  obtained  by  integrating 
acceleration  twice,  then  the  maximum  permissible  constant  error  in 

acceleration  is  2f /At^.  Using  a  value  of  100  ft  for  i  (as  specified  in 
in  accuracy  requirements  for  some  tests),  permits  computation  of 
numerical  values  for  the  maximum  constant  velocity  and  acceleration 
errors  as  functions  of  the  integration  time  At.  Such  relations  are 
shown  by  Figures  1  and  2. 

Constant-Relative -Velocity  Assumption 

A  relative  trajectory  has  a  greatly  simplified  analytic  form  when 
the  assumption  of  linear  relative  trajectory  (or  constant  relative 
velocity)  is  justified.  This  assumption  will  be  justified  if  the  dis¬ 
placement  due  to  relative  acceleration  is  sufficiently  small  to  be 
ignored  for  the  time  interval  of  interest. 


Integration  Time  At,  sec 


Figure  I.  Constant  Velocity  Error  for  100-ft  Accumulated 
Position  Error  vs  Integration  Time  At 


Acceleration  Error  for 


Integration  Time  At,  sec 


Figure  Z.  Constant  Acceleration  Error  for  100-ft 
Accumulated  Position  Error  vs  Integration  Time  At 
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A  typical  scoring  requirement  specifies  that  an  accuracy  of  ±  £ 
ft  in  position  is  desired  for  all  ranges  R  :<  Maximum  and 

minimum  closing  rates  are  also  given,  and  the  time  interval  over 

which  scoring  is  desired  is  then  bounded  by  2R,  ,  /B,.  S  At  s 

,  ^  ’  Max  Max 

2R,,  /R,,.  . 

Max  Min 

For  the  problem  of  scoring  an  attack  on  an  air  vehicle  at  any 

altitude  between  70,  000  and  200,000  ft,  the  requirements  are,  £  =  100ft, 

,  =  3000  ft,  and  6000  <R  <  10,000  ft/sec.  The  smallest  time 

Max 

interval  of  interest  for  this  type  of  attack  is  thus  6000/10,  000  =  0.6  sec; 
the  largest  time  interval  would  be  6000/6000  =  1.0  sec. 

The  requirements  for  scoring  attacks  on  ICBM's  are 

2000  ft  and  £  =  iOO  ft.  A  study  of  the  possible  tactics  for  such  attacks 
indicates  that  they  are  likely  to  be  nose  attacks,  and  the  closing  rates 
will  thus  be  large.  If  the  smallest  closing  rate  is  taken  to  be 
10,  000  ft/sec,  the  maximum  time  interval  for  scoring  will  be 
4000/10, 000  =0.4  sec. 

The  last  problem  to  be  considered  is  that  of  determining  the 

trajectory  of  a  missile  fired  from  one  satellite  against  a  second 

satellite.  It  is  specified  that  R.  .  =  500  ft,  and  £  should  be  about 

^  Max 

15  ft  if  the  accuracy  is  to  be  comparable  to  that  required  for  the  air- 
vehicle  problem.  Closing  rates  would  be  enormous  for  head-on 
attacks  in  this  case  (of  the  order  of  50,  000  ft/sec),  and  it  seems  likely 
that  the  much  lower  rates  associated  with  tail  attacks  would  be  desir¬ 
able  for  accuracy.  A  value  of  2000  ft/sec  is  not  unreasonable  for  the 
minimum  closing  rate;  this  corresponds  to  a  maximum  time  interval 
of  1000/2000  =  0.  5  sec. 

The  maximum  error  in  position  due  to  the  assumption  of 
constant  velocity  is  the  second  integral  of  the  relative  acceleration 
over  the  time  interval.  If  a  maximum  of  Ng's  relative  acceleration  is 
effective  during  the  attack,  the  resulting  range  error  will  be  no  more 

than  AR  =  0.  5  Ng  (At)^. 

If  the  maximum  position  error  is  100  ft, 

iOO  >  16.  1  N(At)^  or  N  <6.211/(At)^ 

which  relates  the  permissible  average  or  constant  relative  acceleration 
to  the  time  interval  of  the  attack  when  the  maximum  error  in  position 
(due  to  the  linear  relative  trajectory  assumption)  is  taken  to  be  100  ft. 

A  graph  of  this  relation  is  shown  in  Figure  3.  If  the  permitted  error  £ 
is  not  iOO  ft,  values  from  this  curve  should  be  multiplied  by  (£/ 100). 
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Since  the  time  interval  to  be  considered  is  related  to  R,  ,  and 
.  ^  Max 

R,  N  can  be  related  to  R  instead  of  At.  For  100-ft  allowable  error 

and  for  R,,  =  3000  ft,  the  relation  is 

Max 

100  >  16.  1  N(6000/R)^  or  N  <R^/5. 796  X  10^ 

This  relation  is  plotted  in  Figure  4.  If  E  is  not  100  ft,  values  from 
this  curve  can  be  corrected  as  before  by  the  factor  (E/100). 

It  will  be  noted  that  R  is  smallest  and  At  is  largest  for  tail 
attacks.  These  attacks  take  longer  times  to  be  completed,  other 
factors  being  equal,  so  that  relative  acceleration  acts  over  a  longer 
time  interval  to  cause  displacement.  Thus  the  permissible  relative 
accelerations  are  smallest  for  tail  attacks.  This  does  not  deny  the 
validity  of  the  assumption  for  certain  tail  chases,  however.  In 
particular,  it  is  likely  to  be  valid  for  the  satellite-vs  -  satellite  problem. 

An  average  relative  acceleration  of  20  g's  can  cause  a  displace¬ 
ment  of  100  ft  in  0.5573  sec.  If  the  range  interval  of  interest  is 
R  <  2000  ft,  then  the  linear-relative-trajectory  assumption  is  valid  for 
closing  rates  greater  than  or  equal  to  4000/0.5573  =  7177  ft/sec.  The 
assumption  of  a  linear  relative  trajectory  over  these  ranges  is  thus 
apparently  valid  for  all  nose  attacks  on  satellites  and  ICBM's. 

If  the  target  and  weapon  have  constant  velocities  relative  to  a 
third  observation  station,  they  have  constant  velocities  relative  to  each 
other.  The  converse  is  not  necessarily  true;  they  may  have  constant 
velocities  relative  to  each  other,  to  give  a  linear  relative  trajectory, 
while  at  the  same  time  both  target  and  weapon  may  have  equal  non¬ 
zero  accelerations  with  respect  to  the  third  observation  point.  In  this 
latter  case,  the  trajectories  of  target  and  weapon  with  respect  to  this 
observation  station  are  not  linear. 

Relative  Motion  of  Munition  V»'ith  Respect  to  Target  for  Constant 
Relative  Velocity 

If  the  velocity  of  the  weapon  relative  to  the  target  can  be  well 
approximated  as  constant  in  magnitude  and  direction  during  an  attack, 
the  relative  trajectory  can  be  given  as  a  linear  function  of  time.  If 

R  is  the  range  vector  of  the  weapon  with  respect  to  the  target,  R^  is 

its  derivative  in  any  inertial  (non-rotating)  coordinate  system,  and  R_^ 

is  the  initial  value  of  R,  then 

R  (t)  =  R  +"R  t 
'  o  o 


(1) 
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R,  Thousands  of  Feet  per  sec 


Figure  4.  Relative  Acceleration  for  100-ft  Accumulated  Position 
Error  vs  Closing  Rate,  with  Maximum  Range  3000  ft 
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The  relative  coordinates  in  any  inertial  system  are  then  given  by 


X(t)  =  X  +  X  t 
'  o  o 

Y(t)  =  Y  +  Y  t 
'  o  o 

Z(t)  =  Z  +  Z  t 
'  o  o 

where 


X 


^oW  ■  ^oT 
^oW  -  ^oT' 


(2) 


(3) 


The  subscripts  W  and  T  denote  weapon  and  target  values,  respectively, 
and  subscript  o  denotes  the  initial  point.  The  length  of  the  range 
vector  is  given  by 


R  =  (X^  +  Y^  +  Z^)" 


R 


(X^  +  Y^  +  Z^)  +  2t{X  X  +Y  Y  +Z  Z  )  +  t^d^  +  Y^  +  Z^)  i  ^ 

r\  r\'  '  r\  r\  r\  r\  r>rs'  'O  O  0^1 


O  O  O  O  0  0 


R  =  (R  •  R)^  =  (R^  +  2R  •  R  t  -1  R^  t^)^ 
'  '  'o  OO  0' 


(4) 


If  R^  is  constant  and  is  aligned  along  R^,  R  is  a  linear  function  of 

the  time,  and  a  ^irect  hit  would  result  for  positive  closing  rates;  also 
the  direction  of  R  remains  constant  during  the  approach  of  the  weapon 

to  the  target.  On  the  other  hand,  if  R^  is  not  aligned  with  R^,  the 

magnitude  of  R  is  not  a  linear  function  of  the  time,  and 


dR 

~5r 


R  •  R  +  R^  t 
o  o  o 

R 


(5) 


The  quantity  dR/dt  is  of  maximum  magnitude  at  extreme  range 
and  decreases  with  time  to  the  minimum  range  point,  after  which  it 
again  increases  in  magnitude.  In  all  cases,  dR/dt  =  0  at  the 
minimum  range  point.  This  gives  a  convenient  method  for  finding 
the  minimum  range  point.  Since  the  numerator  of  equation  (5)  is 
zero  at  minimum  range. 
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=  -R  .  R  /R^ 


R,  ..  =  R  -  R  (R  •  R  )/R^ 

Min  o  o  '  o  0  0 


R 


Min 


R"^  -  (R  .  R  )^/R'^ 
o  '  o  o  o 


(6) 


Observe  also  that  if  R^  is  not  aligned  along  R^,  the  relative 

range  R  is  no  longer  constant  in  direction.  The  minimum  angular 

rate  of  rotation  of  R  then  occurs  at  maximum  range,  and  the  maximum 

angular  rate  of  rotation  of  R  occurs  at  the  minimum  range.  Since  R 

is  the  sum  of  R  and  R  t,  the  relative  motion  occurs  in  the  plane  of 
.  o  o  .  ^ 

R^  and  R^.  In  this  plane  the  relative  velocity  R  has  components 

along  R  and  normal  to  it,  with  the  sums  of  the  squares  of  these 

■**  — 

components  being  equal  to  the  magnitude  of  (R  •  R).  At  the  point  of 
closest  approach  (R  >  =  0  so  that  the  relative  velocity  is  normal 

to  The  velocity  normal  to  R  may  in  general  be  represented  as 

the  vector  cross  product  of  the  rotational  rate  vector  of  R  and  R  itself; 
i.  e.  ,  the  component  of  R  normal  to  R  is  given  by  X  R,,  where 
is  the  vector  angular  rotation  rate  of  R  and  is  normal  to  both  Rand 
R.  Actually, 


R  X  R 


R 


(7) 


This  angular  rate  has  the  maximum  value 


R 


o 


'R 


Max 


R 


Min 


.  1 

R^  -  (R  •  R 

O  '  O  O  '- 


(8) 


This  maximum  angular  rate  tends  to  <»  as  tends  to  zero,  so 

that  large  angular  rates  are.  to  be  expected  for  near  misses  at  the  point 
of  closest  approach.  The  angular  rate  of  R,  given  by  ,  may  be 
resolved  into  components  in  any  convenient  reference  frame,  and 
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azimuth  and  elevation  rates  can  then  be  determined.  Note  that  any 

component  of  must  have  absolute  value  <  .  These  angular 

rates  are  those  of  the  weapon  as  viewed  from  the  target  and  expressed 
in  an  inertial  reference  frame.  The  angular  rates  of  the  weapon  and 
of  the  target  as  seen  from  a  third  reference  point  may  be  considerably 
different  from  these,  but  they  would  be  given  by  similar  relations. 

For  such  linear  relative  trajectories  the  direction  of  R  rotates 

through  almost  180°  from  the  beginning  of  attack  to  final  departure  of 
the  weapon  from  the  immediate  neighborhood  of  the  target.  If  the 
weapon  is  to  be  observed  from  the  target  throughout  the  trajectory,  an 

angular  coverage  of  over  1.80°  would  be  required;  and  if  the  side  on 
which  the  weapon  passes  the  target  is  not  specified  in  advance, 

360°  coverage  should  be  provided.  However,  if  viewing  the  target  up 
to  the  point  of  closest  approach  is  sifficient,  only  half  this  angular 
coverage  would  be  required.  Extrapolation  or  integration  could  be 
used  to  compute  portions  of  the  trajectory  outside  the  field  of  view. 

The  propagation  of  small  errors  in  these  vector  and  analytic 
relations  may  be  studied  by  taking  the  differentials.  Thus, 

dR  =  dR  +  tdi  +  R  dt  (9) 

o  o  o  ' 

This  equation  gives  the  effects  of  small  errors  made  in  measuring 
R^,  R^,  and  t  upon  the  resulting  error  in  R.  Normally,  dt  will  be 
negligibly  small  in  comparison  to  other  errors  in  scoring  tests.  Taking 
the  components  of  dR  gives 


dX  =  d  1  (X^.^  -  X^.j,)  +  t(X^.^  -  X^.j,) 


dX  =  dX^^-dX^^  +  (X^^-X^^)dt  +  t(dX^^-dX^^) 

with  similar  expressions  for  dY  and  dZ. 

The  error  in  the  measured  value  of  R  is  given  by 
dR  =  (Xdx  +  Ydy  +  ZdZ)/R 


(10) 


(11) 
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The  error  in  computed  from  equation  (6)  is  given 

approximately  by 

dR...  =  dR  -  d  |r  (R  •  R  )/R  .  R 


R  ( 

o' 


o'/%  •  ^oj 


dR^,,  =  dR 

Min  o 


dRp  (Rq  °  Rq)  °  ptdRo  •  Rq)  +  R^(R^  •  dR^) 

R^ 


Rp  (^o-  V(2dR^.  R^) 

•  4 

R 

o 


(12) 


If  the  angular  rotation  rate  is  computed  from  equation  (7), 
the  differential  error  is  given  by 


di7 


R 


dR  X  R  +  R  X  dR  2(R  X  R)(R  •  dR) 


R 


R 


(13) 


Velocity  and  Acceleration  of  a  Rocket  (Atmospheric  Drag  and  Gravity 
Effects  Neglected) 


Simple  relations  may  be  obtained  to  approximate  the  velocity  and 
acceleration  functions  of  a  rocket  when  atmospheric  drag  and  gravity 
effects  are  neglected.  These  relations  will  be  valid  when  the  thrust  is 
large  compared  to  the  neglected  effects.  It  will  be  assumed  that  the 
thrust  is  aligned  with  the  velocity  vector. 

The  general  equation  for  the  acceleration  reduces  in  this  case  to  . 


mV  =  -Ng  m  =  thrust 


(14) 


where  the  dot  denotes  the  time  derivative,  and 
V  =  velocity,  ft/sec 

m  =  mass  of  rocket  and  remaining  fuel,  slugs 
m  =  negative  of  fuel  burning  rate,  slugs/sec 

N  =  specific  impulse,  Ib-sec/lb  or  sec 

2 

g  =  gravitational  constant,  ft/sec 
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This  equation  may  be  integrated  to  give  the  velocity  increment. 


dV 


-Ng 


dm 

m 


AV  =  V-V  =  Ng  In  -  ^ 

°  1  -  (m^/m^) 

where  m^  is  the  mass  of  fuel  expended,  and  m^  is  the  initial  total  mass 
of  rocket  and  fuel. 

Figure  5  is  a  plot  of  AV  vs  m^/m^  for  three  values  of  N.  Since 

equation  (15)  shows  that  the  velocity  increment  is  linear  in  N,  the 
value  of  AV  for  any  N  is  {N/300)AV(300,  m^/m^}  where  AV{300,  m^/m^) 

is  the  value  of  AV  read  from  the  curve  for  N  =  3  00.  Figure  6  is  a  plot 
of  AV/Ng  vs  m^./m^,  which  is  a  dimensionless  curve  serving  for  all 

values  of  the  variables. 

If  the  burning  rate  m  is  constant,  m  =  m^  +  mt,  and  equation  (14) 
becomes 


-Ng  m 

m  +  mt 
o 


which  can  be  written  in  the  dimensionless  form 


Vm 

o 


1 


Ngm  1  +  {rn/m^)t 


(16) 


^  9  # 

Figure  7  shows  a  plot  of  (-Vm^/Ngm)  vs  (-mt/m^).  If  N,  ni^,  and 

m  are  given,  acceleration  V  is  a  function  of  motor  burning  time  can  be 
obtained  from  this  curve  for  any  time  0  <  -mt/m^Sm^/m^  where  m^ 

is  the  total  mass  of  the  fuel. 


The  velocity  increment  as  a  function  of  motor  burning  time  for 
constant  fuel  burning  rate  may  be  similarly  obtained.  In  this  case 
equation  (15)  becomes 


AV  =  V  -  V^  =  Ng 


In 


1 


1  + 


mt 

m 

o 


(17) 
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Converting  this  equation  to  dimensionless  form  gives 


(18) 


Figure  8  is  a  plot  of  this  relation,  from  which  the  velocity  incre¬ 
ment  as  a  function  of  the  burning  time  may  be  found.  This  curve  is  the 

same  as  Figure  6  with  -mt  substituted  for  m^. 

The  displacement  or  change  in  position  due  to  motor  burning  at 
a  constant  rate  may  be  obtained  by  integrating  equation  (17)  with 
respect  to  time.  If  S  is  the  total  displacement  during  motor  burning, 
then  the  displacement  due  to  acceleration  is  S  -  V  t.  Integration  of 
equation  (17)  gives 


-m  (S  -  V^t) 


\  my  \  m  y  m 

'  r\^  '  r\  '  r\ 


(19) 


Figure  9  is  a  plot  of  this  dimensionless  displacement  due  to  motor 
burning  vs  the  dimensionless  time  (-mt/m^). 

Figure  10  is  a  dimensional  plot  of  the  displacement  vs  time  for 

a  hypothetical  rocket  with  constant  burning  rate  m  =  -m  /ISO  slugs/sec 
and  specific  impulse  N  =  300  sec.  The  curve  ends  at  t  =  135  sec 
corresponding  to  a  fuel-mass  ratio  of  0.  9.  The  acceleration  increases 
with  time  to  a  maximum  of  20  g's  at  burnout.  It  appears  that  a 
constant  burning  rate  is  not  the  most  efficient  way  to  program  motor 
burning  because  larger  displacements  and  the  same  final  velocity 
result  if  the  larger  accelerations  occur  at  the  beginning  of  the  flight. 

Relative  Trajectory  with  Non-Zero  Relative  Acceleration 

The  velocity  and  position  of  a  weapon  relative  to  a  target  may  be 
obtained  by  integrating  the  relative  acceleration.  If  the  relative 
acceleration  is  n(t),  then  the  relative  velocity  is  given  by 

V(t)  =  +  y^n(ti)  dt^  (20) 


The  relative  displacement  of  the  weapon  with  respect  to  the  target  is 
then  given  by 


Lmensionless  Motor  Burning  Time,  (-mt/m^ 


.  Dimensionless  Velocity  Increment  as  a  Fun 
inless  Motor  Burning  Time  for  Constant  Burn: 
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‘?uatu0D-ei;dsta 
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S(0  .  t  Y(t)  d,  =  +  V„t  t  dt^  „(.,)  d., 


(21) 


The  relative  velocity  and  displacement  cannot  generally  be  given 
as  closed-form  functions  of  time;  but  if  the  relative  acceleration,  n(t), 
is  a  linear  function  of  time,  the  integrations  can  be  performed.  Let 
n(t)  =  A  +  Bt;  then 


V(t) 


=  V  +  r  (A  +  Bt,)dt,  =V  +  At  +  Bt^/2 
o  Jq  1'  1  o 


S(t)  =  +  V^t  +  AtV2  +  BtV^ 


(22) 


Higher-order  terms  could  be  added  to  these  equations,  but  it 
seems  likely  that  the  scoring  encounters  with  which  this  study  is 
concerned  will  be  characterized  by  accelerations  which  are  no  worse 
than  linear  functions  of  time.  If  it  can  be  assumed  that  the  underlying 
function  represented  by  the  acceleration  data  is  a  polynomial  in  time 
of  degree  n  or  less,  where  n.  is  known,  the  curve-fitting  coefficients 
A,  B,  etc.  do  not  have  to  be  determined.  Instead,  smoothing  and 
interpolation  formulas  can  b«  developed  which  will  perform  the 
integration  and  give  the  trajectory  directly.  Reference  2  gives  such 
formulas  for  polynomials  of  degree  3  or  less. 

Equations  (22)  represent  the  velocity  and  displacement  for  an 
acceleration  which  is  a  scalar  function  of  the  time.  Similar  relations 
would  be  obtained  for  each  of  the  rectangular  coordinate  components 
if  a  vector  acceleration  were  assumed. 

Reference  i.  "The  Trajectory  of  a  Rocket  with  Thrust,"  R.  E. 
Struble,  C.  E.  Stewart,  and  J.  Granton,  Jr.,  Jet  Propulsion,  Vol.  28 
No.  7,  July  1958,  pp.  472-478. 

Referense  2.  Mathematical  Processes:  Sets  of  Coefficients  for 
Data  Processing  Formulas,  Evelyn  Welborn,  MPRL  Report  475, 
Military  Physics  Research  Laboratory,  The  University  of  Texas, 
Austin,  Texas,  5  May  1959. 
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APPENDIX  2 

EFFECTS  OF  OWN-SHIP  ANGULAR  MOTION  ON  RELATIVE 

TRAJECTORIES 


Introduction 

The  desired,  objective  in  taking  data  on  missile-target  scoring 
encounters  is  to  secure  enough  information  to  establish  within 
permissible  error  limits  the  trajectory  of  the  weapon  in  the  vicinity  of 
and  with  respect  to  the  target.  It  is  highly  desirable  that  this  relative 
trajectory  be  referred  to  a  stable  coordinate  system  whose  approximate 
absolute  orientation  is  known.  If  the  trajectory  scoring  method  is 
based  on  assumptions  about  the  nature  of  the  relative  trajectory,  the 
data  must  be  measured  in  (or  at  least  be  converted  to)  a  stable 
system.  A  stable  coordinate  system  is  here  defined  as  one  which  does 
not  rotate  appreciably  during  the  scoring  encounter. 

It  will  be  necessary  in  many  instances  to  obtain  trajectory  data 
by  use  of  instruments  rigidly  attached  to  the  target  or  missile;  the 
data  taken  will  be  referred  to  body  axes  of  the  parent  object  and  will 
thus  be  taken  with  respect  to  a  rotating  coordinate  system.  The  own- 
ship  angular  motions  would  then  introduce  errors  into  the  data  and 
into  the  relative  trajectory  obtained  from  the  data;  these  errors  must 
be  taken  into  account  and  corrected  for  when  they  are  sufficiently 
large. 


The  linear  velocities  and  accelerations  of  the  own-ship  center  of 
mass  affect  the  relative  trajectory,  of  course,  but  they  do  not  of 
themselves  introduce  errors  into  observations  taken  from  one  of  the 
objects.  Errors  will  be  introduced  if  instruments  which  are  not 
intended  to  measure  these  quantities  are  sensitive  to  them  and  if  the 
scoring  method  is  based  on  invalid  assumptions  about  the  nature  of  the 
relative  trajectory.  These  problems  should  be  considered  in  the  design 
of  the  scoring  system  but  are  not  of  present  concern. 

The  purposes  of  this  paper  are  to  examine  the  effects  of  own- 
ship  angular  motions  on  the  data,  to  establish  criteria  for  neglecting 
these  motions,  and  to  discuss  methods  of  taking  them  into  account 
when  they  may  not  be  neglected. 

Trajectory  Errors  Due  to  Own-Ship  Angular  Motions 

Figure  1  shows  a  stable  rectangular  Cartesian  coordinate 
system  X^,  Y^,  with  origin  at  the  center  of  mass  of  the 

instrumented  object;  the  rotating  system  X,Y,Z  is  fixed  to  the  instru¬ 
mented  object  and  is  shown  with  its  origin  O  coinciding  with  O^.  The 

vector  range  from  O  to  the  observed  object  at  point  P  is  denoted  by  R. 
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The  azimuth  and  elevation  angles  of  R,  measured  with  respect  to  the 
rotating  system,  are  denoted  by  A  and  E,  respectively.  Angle  A  is 
measured  in  the  X-Z  plane  from  the  X  axis  to  the  plane  o^R  and  Y  and 
is  positive  as  sho_wn;  angle  E  is  measured  in  the  plane  of  R  and  Y  from 
the  X-Z  plane  to  R  and  is  positive  as  shown. 

Figure  2  shows  the  situation  when  the  origins  are  separated; 

here  D  is  the  vector  from  O  to  O,  and  R  is  the  vector  from  O  to  P. 

o  o  o 

Since  measurements  may  be  taken  at  two  or  more  stations,  it  may  be 

necessary  in  practice  to  consider  two  or  more  rotating  systems; 

(X^,Yj^,Z^),  (X^,,  Y^.  Z.^),  etc.  with  origins  at  Oj^,  O^,  etc.  If  there  is 

no  b_ending,  the  rotating  systems  may  all  be  parallel,  but  the  vectors 
Dp  ,  etc.  will  all  be  different.  For  simplicity,  this  discussion  will 

be  limited  to  a  single  rotating  system;  the  extension  to  multiple  systems 
is  obvious . 

Bending  can  occur  as  a  result  of  aerodynamic  loading,  warping 
of  the  structure  due  to  heating  or  cooling,  etc.  ;  it  will  have  to  be 
measured  separately  from  the  own-ship  angular  motions.  The  effects 
of  bending  will  be  to  move  the  origin  O  and  twist  the  X,Y,Z  system. 

The  coordinate  transformation  to  correct  for  bending  will  therefore 
contain  translation  as  well  as  rotation  terms  and  could  be  developed 
from  a  diagram  similar  to  that  of  Figure  2.  The  dynamic  effects  of 
bending  may  be  small  enough  to  neglect;  but  if  the  angular  velocity  and 
angular  acceleration  of  the  coordinate  system  are  measured  at  or 
near  origin  O,  no  assumption  need  be  made  about  this  point.  The 
change  in  the  length  of  D  is  probably  negligible,  but  it  would  have  to  be 
measured  or  determined  from  a  calibration  if  it  is  not. 

Values  of  the  coordinates  of  point  P  in  the  X^,  Y^,  Z^  system 

will  be  determined  by  applying  in  succession  the  coordinate  transfor¬ 
mations  for  bending  and  for  own- ship  angular  motions  to  the  coordinates 
measured  in  the  X,Y,Z  system;  but  if  the  method  of  scoring  assumes 
no  bending,  some  preliminary  calculations  will  be  required.  An 
example  of  this  is  the  two- station  angle -only  triangulation  system  in 
which  it  is  assumed  that  O  is  a  fixed  (known)  distance  from  the  second 
station  and  that  X,  Y,  Z  has  a  specified  orientation  with  respect  to  the 
line  connecting  the  stations;  similar  assumptions  are  made  about  the 
coordinate  system  at  the  second  station.  If  significant  bending  occurs, 
it  will  be  necessary  to  determine  new  standard  orientations  for  the 
stations  and  then  to  correct  the  angular  data  and  the  length  of  the  base 
line.  After  the  corrected  data  have  been  used  to  compute  the  range, 
the  coordinate  transformations  for  bending  and  for  own- ship  angular 
motions  would  be  applied,  one  after  the  other,  to  the  values  of  X,  Y,  Z 
determined  from  the  computed,  range  R  and  the  original  uncorrected 
angles . 
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The  angles  used  to  specify  the  orientation  of  X,  Y,Z  with 
respect  to  X  ,  Y^,  may  vary  from  one  application  to  the  next,  and 

it  is  therefore  impossible  to  write  specific  coordinate  transformations 
for  general  use.  It  can  be  said  that  if  a  stabilized  platform  is  used, 
the  orientation  angles  will  be  defined  by  the  gimbal  system;  similarly, 
the  bending  angles  are  determined  by  the  measurement  system. 

The  coordinate  transformation  from  the  X,  Y,  Z  system  to  the 

X  ,  Y  ,  Z  system  can  be  written,  in  the  form 
o  o  o 

^o  "  ^^X  X^'''^X  Y^'^'^X  Z^^  ■  “^X  X^  Ej^cos  Ap 
o  o  o  o 

Yo  =  (Cy  x^+'^Y  y’^+  ‘^Y  Z^'>  -  Cy  y  ° 

o  o  o  o 

Z  =  (C^  Y^  *^Z  Z^^  ’  *^Z  Z  ^  "^D 

o  “o  o  o 

where  Cy^  y,  etc.  are  direction  cosines.  Since  the  components 

o  o 

of  D  in  the  X,  Y,Z  system  are  constant  (neglecting  be_nding  which  can  be 

treated  separately),  it  is  more  convenient  to  expand  D  in  this  system 

than  in  the  X  ,  Y  ,  Z  system.  The  general  form  of  the  coordinate 
o  o  o  •  ° 

transformation  for  bending  will  be  identical  to  equations  (1)  except  for 

changes  in  symbols.  When  these  transformations  are  applied,  the 

translation  and  rotation  terms  c_an  be  evaluated  separately  so  that 

corrected  components  of  D  and  R  in  the  X  ,  Y  ,  Z  system  will  be 
,  .  .  1  o  o  o 

obtained. 

The  true  velocity  and  acceleration  vectors  R  =  dR/dt  and 

R  =  d^R/dt^  transform  by  equations  (1)  with  the  translation  terms  set 
equal  to  zero  (by  the  product  of  two  such  transformations  if  bending 
occurs),  but  it  should  be  noted  that  these  vectors  are  not  the  ones 
which  would  be  measured  by  instruments  fixed  in  the  X,  Y,  Z  system. 
This  can  be  shown  by  differentiating  the  equation 

R  =  D  +  R 
o 

It  is  convenient  to  consider  the  derivatives  of  D  and  R  separately. 
Expressions  for  the  derivatives  of  R  will  be  developed  first. 

If  i  ,  j  ,  and  k  ar_e  unit  vectors  in  the  X,  Y,  Z  system,  the 
relative  range  vector  R  will  be  denoted  by 

=  Tx  +Ty  +  kZ 


R 


(2) 
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Differentiation  gives 

R  =  (Tx  +Ty  +  kZ)  +  (fx  +Ty  +  kz)  (3) 

•  0  »  * 

R=V+wXR=R-D=V-D  (4) 

o  o 

Z.  ' 

since  it  can  be  shown  that  i  =  ti3Xi,  j  =  coXj,  and  k  =  to  X  k  where 

CO  =  i  coj^  +  j  cOy  +  kco^  (5) 

is  the  angular  velocity  of  system  X,  Y,  Z  with  respect  to  system 

X  ,Y  ,Z  ,  Note  that  the  equation 
o  o  o 


R 


dR 

dt 


CO  X  R 


_a_ 

dt 


+  CO  X 


R 


(6) 


expresses  a  general  operator  relation  applying  to  any  vector  measured 
v/ith  respect  to  the  X,Y,Z  system.  Also  note  that  the  assumption  that 
X^,  Y^,Z^  is  a  stable  system  is  unnecessary  to  this  developm.ent;  the 

equations  apply  to  any  two  systems  which  have  relative  rotation. 

The  subscript  "o"  will  be  used  to  denote  quantities  measured 
with  respect  to  the  X^Y^,  Z^  system;  thus  V  is  the  velocity  of  point  P 

as  observed  in  the  X,Y,Z  system,  and  =  R^  is  the  velocity  of  P 

as  observed  in  the  X  ,Y  ,Z  system.  For  convenience,  a  dot  over 

the  symbol  has  been  employed  instead  of  the  operator  d/ dt  to  denote 
the  total  time  derivative  (or  the  time  rate  of  change  of  the  quantity 
■—v/ith  r-eepect  to  the  desired  reference  system  X^,  Y^.,,  Z^);  the  operator 

8/9t  gives  the  time  rate  of  change  measured  with  respect  to  the 
X,  Y  ,  Z  system.. 

Equation  (4)  shows  that  an  observer  in  the  X,  Y,  Z  coordinate 
system  who  did  not  all^w  for  own^-sh^  rotation  would  report  the 
velocity  of  pointJP  as  V  =  i^l  +  j  Yj^+  1^4,_which  is  smaller  than  the 
desired  vector  V  by  the  amount  D  +  coX  R.  If  the  same  observer  were 

to  determine  the  range  vector  at  any  time  t  by  integrating  the  velocity 
compo^nts  Y,  and with  the  required  initial  conditions,  he  would 
obtain  R=iX  +  jY  +  kZ,  which  is  the  correct  value  for_his  rotating 
coordinate  system  but  is  not  equal  to  the  desired  vector  R^. 

The  accumulated  error  in  position  would  be  given  by  the  time 
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integral  of  (D  +  w  X  R).  If  the  integral  of  D  is  taken  to  be  negligible, 

the  position  error  due  to  own-ship  rotation  would  be  bounded  by 

I  AR  I  I  CO,  ,  R,  ,  At  I,  where  At  is  the  time  interval  involved.  This 

bound  for  position  error  due  to  rotation  of  the  reference  coordinate 
system  applies  whether  the  relative  trajectory  is  obtained  by  direct 
measurement  or  by  integrating  the  relative  velocity  components.  This 
bound  would  ordinarily  be  used  to  determine  the  error  due  to  own- ship 
angular  motion,  but  it  could  also  be  used  to  determine  the  error  due 
to  using  a  stabilization  scheme  which  allows  some  rotation  of  the  ,^,xes 
with  respect  to  inertial  space. 


The  effects  of  angular  rotation  of  the  reference  coordinate 
system  upon  relative  acceleration  will  now  be  examined.  Differentiation 
of  equation  (3)  gives 


R  =  (iX  +  jY  +  kZ)  +  2{iX  +  j  Y  +  kZ)  +  (iX  +  J  Y  +  kZ) 


R=a+2coXV+(iX+jY  + 


•  > 

kZ) 


(7) 


where  a  is  the  relative  acceleration  of  point  P  as  measured  in  the. 

X,  Y,Z  system  regarded  as  non- rotating.  If  equation  (4)  is  differen¬ 
tiated,  use  being  made  of  equation  (6),  there  results 

*•  « 

R=a+2wXV+coXR  +  coX(a)XR)  (8) 

If  D  is  zero  or  small  compared  to  R,  the  last  three  term.s  in  the  right 
member  of  equation  (8)  represent  the  error  that  is  made  in  determining 
the  relative  acceleration  when  the  own-ship  rotation  is  neglected.  It 
may  also  be  shown  by  differentiating  equation  (5)  that  the  angular 
acceleration  of  the  coordinate  system  as  measured  in  the  rotating 

coordinate  system  is  the  same  vector  as  that  measured  in  the  X  ,  Y  , 

o  o’ 

Z  system: 
o 


d(o  -f  •  ,  -f  .-y-*  _8co 

■gr  “  ’-“x  ^  '^“y  ^“z  "  irr 


(9) 


The  effect  upon  the  resulting  trajectory  of  the  error  in  relative 
acceleration  due  to  the  rotation  of  the  reference  coordinate  system 
may  be  obtained  without  using  equation  (8).  If  the  incorrect  acceleration, 

.  -  r,** 

a  =  iX  +  jY  +  kZ,  is  integrated  with  respect  to  the  time,_^us_e  being 
made  of  the  proper  initial  conditions,  the  result  will  be  V  =  i!^  +  jY  +  kZ . 
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The  error  made  in  measuring  relative  acceleration  thus  leads  to  the 
previously  given  errors  in  velocity  and  position. 

Equations  for  the  derivatives  of  D  will  now  be  developed,  if 
where  1^^  is  a  unit  vector  in  the  direction  of  D,  it  will  be 
seen  from  the  previous  discussion  that 

4 

D  =  X  Tjj)  =  +  Wp  X  D 

u_*.. 

D  =  l^D  +  2D(t0j^  X  Ij^)  +  X  D  +  Wq  X  X  D  ) 

where  co  is  the  angular  velocity  of  1  with  respect  to  the  X  .Y  Z 

^  o  o  o 

system,  if  there  is  no  bending,  D  =  D  =  0  and  can  be  determined 
from  the  equation 

where  flp  =  Ip  Op  is  the  angular  velocity  of  the  stabilized  platform  with 

respect  to  the  vehicle.  If  there  is  bending,  D  and  D  are  not  zero  and 
are  not  necessarily  negligible;  moreover,  will  either  have  to  be 

measured  directly  or  the  value  given  by  the  above  equation  will  have  to 
be  corrected  by  use  of  data  from  the  instruments  which  measure  the 
bending.  Numerical  differentiation  will  be  needed  if  the  instruments 
do  not  measure  the  correction  and  the  derivatives  directly. 

The  equations  for  the  velocity  and  acceleration  of  point  P  as 
observed  in  the  X^,  Y^,  Z^  system  can  now  be  written  by  use  of 

equations  (4)  and  (8): 

^o  =  \  =  D  +  iiTxR  +  V  (10) 

-*•  M.  ^  ^  ^  ^ 

^o  ”  ^o  ”  D+coX(wXR)  +  wXR  +  2uXV+a 

A 

where  D  and  D  are  given  by  the  equations  of  the  previous  paragraph. 

If  point  P  is  motionless  with  respect  to  the  X,  Y,  Z  system,  V  and  a' 
are  zero.  The  remaining  terms  in  the  right  members  of  equations  (10) 
are  said  to  be  the  velocity  of  transport  and  the  acceleration  of  transport 
of  the  X,  Y,  Z  system  with  respect  to  the  X^,  Y^,  Z^  system. 

It  is  assumed  in  the  foregoing  discussion  that  the  vector  quantities 
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V  =  iX  +  jY  +  kZ  and  a=iX  +  jY  +  kZ  are  determined  without 
instrumentation  error  either  by  direct  measurement  or  by  numerical 
differentiation  of  the  functions  X(t),  Y(t),  and  Z(t).  Direct  measure¬ 
ments  could  be  made  by  use  of  electromagnetic  radiation,  for  example. 

The  relative  acceleration  of  one  vehicle  with  respect  to  another 
could  be  determined  indirectly  by  use  of  three  inertial  accelerometers 
in  each  vehicle;  each  set  of  accelerometers  would  be  mounted  on  a 
stabilized  platform  with  their  sensitive  axes  mutually  orthogonal. 
Similarly,  the  change  in  relative  velocity  as  a  function  of  time  could  be 
determined  by  use  of  three  inertial  velocimeters  in  each  vehicle.  In 
either  case,  the  stabilized  platforms  should  have  known  and  preferably 
identical  orientations  with  respect  to  an  external  reference;  also  the 
initial  relative  position  and  velocity  of  one  vehicle  with  respect  to  the 
other  would  have  to  be  determined  by  another  method  before  the  data 
could  be  integrated  (numerically)  to  obtain  the  relative  trajectory.  In 
general  it  would  also  be  necessary  to  know  the  gravitational  attraction 
acting  on  each  vehicle;  but  if  only  short-duration,  short-range  relative 
trajectories  are  of  interest  and  if  the  initial  conditions  for  the  integra¬ 
tion  can  be  obtained  during  the  intercept,  this  requirement  can  be 
waived  without  introducing  significant  errors  into  the  computed  tra¬ 
jectory.  Some  of  these  statements  are  fairly  obvious,  but  others 
require  some  amplification. 

An  inertial  accelerometer  measures  the  force  acting  on  a  known 
mass  suspended  in  such  a  way  that  it  has  one  degree  of  translational 
freedom.  Some  of  these  instruments  do  not  indicate  acceleration  but 
give  a  pulse  for  each  incremental  change  in  velocity  AV  =  at  so  that  a 
count  of  the  pulses  gives  the  change  in  velocity  since  the  beginning  of 
measurement;  these  instruments  are  velocimeters  even  though  the 
absolute  velocity  is  not  determined  unless  the  initial  velocity  happens 
to  be  zero.  In  general  it  will  be  necessary  to  measure  V  (and  R 
as  well)  separately. 

For  determining  how  accelerometer  readings  are  affected  by 
movement  of  the  coordinate  system,  it  will  be  assumed  for  the  moment 
that  X^,Y^,Z^  is  a  true  inertial  system  (motionless  with  respect  to 

the  stars)  so  that  Newton's  second  law  applies  to  accelerations 
measured  with  respect  to  that  system.  The  force  on  a  test  mass  m 
attached  to  the  moving  X.,Y..  ,Z^  system  with  origin  will  then  be 
given  by 

f  ♦  • 

F  =  ma  -  mG,  =  m(D,  +a),  X  (oj,,  X  R,  )  +  u,  X  R,  ) 
o  1  '  1  1  '  1  Im  1  Im 

-  mGj^  =  m(a,j,^  -  G^) 
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where  is  the  angular  velocity  of  the  system  with 

respect  to  inertial  space,  is  the  vector  from  origin  Oj^  to  the 

test  mass,  Dj^  is  the  vector  from.  to  Oj^,  mGj^  is  the  gravitational 
attraction,  and  a^^^  is  the  acceleration  of  transport. 

An  inertial  accelerometer  located  at  the  point  defined  by  R, 

measures  the  component  of  which  is  parallel  to  its  sensitive  axis, 

but  it  also  detects  and  is  biased  by  the  gravitational  field  (Newtonian 
attraction)  at  that  point.  If  the  motion  of  the  accelerometer  (or  of  the 
vehicle  containing  the  accelerometer)  is  unconstrained,_^the  reaction 
of  the  accelerometer  to  the  gravitational  component  of  a.pj  is  equal 

and  opposite  to  the  bias  produced  by  the  gravitational  field.  Thus  the 
inertial  accelerometer  does  not  measure  the  acceleration  due  to  the 
gravitational  field  although  it  can  be  used  to  determine  that  acceleration 
if  the  correct  value  of  a.^.  is  known  or  can  be  computed  or  if  the  motion 

of  the  accelerometer  is  suitably  constrained.  If  G  is  the  gravitational 
field  at  the  point  of  measjarern^nt,  _the  vector  determined  from  acceler¬ 
ometer  readings  will  be  a  =  a_  -  G. 

3.  X 

It  will  be  seen  from  the  second  of  equations  (10)  that  the 
acceleration  of  transport  of  a  second  moving  system  X^,  Y^i  with 

origin  can  be  related  to  the  m.otion  of  the  X^,Yj,Zj^  system  as 

follows: 

• «  * 

-  10^  +  (i)|X  ^12  ^  ^  ^21 

“21  ^  (“21  ^  '^m^  “21  ^  ’^m  “  ®'a2 

where  R^^  vector  from  origin  to  origin  O^,  vector 

from  to  the  point  of  measurement,  is  the  relative  velocity  of  the 

two  systems,  is  the  relative  acceleration,  and  is  the  angular 

velocity  of  the  second  system  with  respect  to  the  first  system.  If  the 
accelerometers  in  the  two  systems  are  mounted  on  stabilized  plat¬ 
forms,  the  vectors  coj^,  w^,  <^2l’  “2I’  ^Im’  ^m  small 

or  zero;  and  terms  containing  Rj^  s-nd  R^j^  will  certainly  be  negligible. 
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Thus 


a2 


al 


=  X  (oij^  X  1^22^  ^ 


“l  ^  ^12  '*' 


^21 


+  a 


21 


(C^-Gi) 


which  is,  except  for  the  last  term,  the  relative  acceleration  of  the 
X^i  Y^i  system  with  respect  to  the  Xj^,  Yj^,  Z  system  corrected  for 

rotation. 

The  quantity  (G_,  -  Gj^)  is  negligible  for  short-duration,  short- 

range  relative  trajectories  but  cannot  be  neglected  when  the  obtainable 
initial  conditions  require  determination  of  the  trajectory  over  extended 
periods  of  time  by  purely  inertial  means.  The  two  coordinate  systems 
are  preferably  parallel  but  obviously  must  be  synchronized  (i.  e.  ,  the 
orientation  of  X^,  Y^i  Z^  with  respect  to  Xj^,Y,^,Z^  must  be  known  for 

all  times  of  interest)  in  order  for  the  subtraction  to  be  performed. 
Each  coordinate  system  could  be  stabilized  with  one  axis  vertical  and 
one  axis  pointing  to  geodetic  north,  an  arrangement  frequently  prefer¬ 
red  for  navigation  or  guidance  purposes,  but  the  two  systems  would 
then  have  a  small  relative  rotation  ^  0).  Assuming  perfect 

stabilization,  the  effects  of  this  angular  velocity  and  of  the  slight  and 
varying  misalignment  of  the  two  systems  would  be  negligible  for  the 
short-range,  short-duration  relative  trajectories  of  present  interest. 

It  would  of  course  be  preferable  in  general  to  have  the  two  coordinate 
systems  non-rotating. 

Criteria  for  Neglecting  Own -Ship  Angular  Motions 


Under  certain  conditions,  own-ship  angular  m.otions  may  be 
small  enough,  or  the  duration  of  the  missile-target  encounter  may  be 
short  enough,  for  these  motions  to  be  neglected.  Criteria  for 
neglecting  own-ship  angular  motions  will  be  developed  in  the  following 
paragraphs  for  the  three  scoring  situations  of  interest  in  this  study. 

It  should  be  remembered  that  these  criteria  can  also  be  used  in  data 
reduction.  Provision  should  be  made  for  correcting  for  own-ship 
angular  motions  unless  it  can  be  demonstrated  that  they  will  always  be 
negligible. 

Three  future  USAF  scoring  requirements  are  specified  by  the 
contract  for  consideration  in  this  study.  The  relative  trajectory  of 
the  missile  with  respect  to  the  target  is  wanted  in  each  case  for  all 
relative  ranges  R  <  The  first  scoring  problem  is  the 

trajectory  of  a  missile  with  closing  rates  of  6000  to  10,  000  ft/sec  at 

altitudes  of  70,  000  to  200,  000  ft  with  R.  ,  =  3000  ft.  The  allowable 

Max 

position  error  is  ±100  ft;  this  is  a  relative  accuracy  of  one  part  in 
thirty. 
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The  error  due  to  angular  velocity  alone  cannot  exceed  the 
total  allowable  error  minus  the  error  from  sources  other  than  angular 
velocity.  The  maximum  position  error  due  to  co  is  given  by  the 
expression  developed  in  the  preceding  section.  The  maximum  time 
duration  of  these  trajectories  is  6000  ft/6000  ft/sec  =  1  sec.  Then  in 
order  for  angular  velocity  to  of  the  reference  system  to  introduce 
negligible  trajectory  error. 


|R.,  w,.  At..  I  =  3000  I  CO,.  I 

'  MaxTvIax  Max'  '  Max' 


100  ft  -  (Position  err_pr 
not  due  to  co 


<  0.  033  rad/sec  -  (position  error  not  due  to  co)/3000 


(li) 


The  fraction  of  the  total  allowable  error  to  be  allocated  to 
angular-velocity  effect  is  arbitrary;  the  choice  will  depend  on  the 
normal  behavior  of  the  scorer  vehmle  and  on  the  size  of  the  remaining 
error.  Thus  the  error  not  due  to  co  may  be  so  large  that  no  allowance 
can  be  made  for  uncorrected  own-ship  motions;  also  the  normal 
angular  rates  of  the  vehicle  may  be  so  larqe  that  there  is  very  little 
probability  of  their  even  being  negligible.  These  remarks  apply  to  any 
of  the  scoring  situations  being  considered. 

The  angular  velocity  criterion  expressed  by  equation  (11)  is 
conservative.  A  statistical  approach  could  probably  be  justified  for 
the  most  part,  but  the  resulting  possibly  more  liberal  criterion  would 
involve  an  arbitrary  relation  between  the  maximum  permissible 
probable  error  due  to  neglecting  co  and  the  error  bound  developed  in 
the  previous  section.  Since  the  criteria  given  here  are  at  the  best 
preliminary  estimates  subject  to  modification  by  experience,  this 
alternate  approach  is  not  attractive. 

For  satellite-vs-satellite  encounters,  scoring  is  required  over 
the  range  interval  of  0  to  500  ft;  the  maximum  permissible  position 
error  is  taken  to  be  one  part  in  thirty  or  about  15  ft.  If  a  minimum 
closing  rate  of  2000  ft/sec  is  assumed,  then  a  maximum  time  interval 
of  1000  ft/2000  ft/sec  =  ^.  5  sec  is  of  interest  for  such  attacks.  Then 
for  the  angular  velocity  co  of  the  reference  coordinate  system  to 
introduce  negligible  trajectory  error, 


R 


Max 


CO.  r  At,,  Is  15  ft  -  (position  error  not  due  to  co) 
Tc/[ax  Max' 


I'^Iax  I  ^  rad/sec  -  (position  error  not  due  to  co  )/250  (12) 


For  scoring  attacks  on  ICBM's  a  range  interval  of  interest  of 
0  to  2000  ft  is  specified.  During  all  three  phases  of  the  ICBM 
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trajectory  (boost,  mid-course,  and  re-entry)  the  attacks  probably 
would  occur  in  the  forward  hemisphere  of  the  ICBM,  and  thus  large 
closing  rates  would  be  expected.  The  smallest  closing  rates  would 
occur  for  attacks  in  which  the  weapon  would  be  deployed  in  the  path  of 
the  on-coming  ICBM  and  would  have  comparatively  little  velocity  of 
its  own.  The  smallest  closing  rate  could  be  expected  to  be  at  least 
of  the  order  of  10,  000  ft/sec,  so  that  a  change  in  range  of  4000  feet 
corresponds  to  an  attack  duration  no  greater  than  0.4  sec_^  If  the 
total  position  error  is  to  be  held  to  ±  100  ft,  the  bound  on  co  would  be 

I^McLxI  ~  rad/sec  -  (position  error  not  due  to  a))/800  (13) 


Elimination  of  Own -Ship  Angular  Velocity  Effects  from  the  Data 

When  the  own-ship  angular  velocity  is  too  large  to  be  neglected, 
as  indicated  by  the  simple  criteria  of  the  last  section,  it  must  in  some 
way  be  taken  into  account.  The  simplest  way,  from  the  viewpoint  of 
data  reduction,  is  to  stabilize  the  reference  axis  system.  This  could  be 
done  by  stabilizing  either  the  entire  vehicle  or  a  platform  carrying 
the  data- gathering  instruments.  Stabilization  of  the  entire  vehicle  would 
not  be  applicable  where  it  interferes  with  the  normal  operation  of  the 
vehicle  in  the  performance  of  its  mission.  It  would  thus  not  be 
applicable  to  scoring  encounters  in  the  atmosphere  or  to  attacks  on  an 
ICBM  but  might  be  applicable  to  satellite-vs  -  satellite  encounters. 
Disturbing  forces  acting  to  change  the  orientation  of  the  satellite  would 
be  4uite  small,  and  stabilization  could  be  maintained  with  little  expendi¬ 
ture  of  fuel  and  power.  Further,  such  stabilization  might  be  needed 
during  attacks  to  minimize  missile-release  problems..  Stabilization  of 
the  instrument  platform  alone  might  be  indicated  for  attacks  on 
instrumented  ICBM's  and  for  attacks  within  the  atmosphere.  Provision 
should  be  made  for  measuring  bending  if  there  is  any  likelihood  that 
it  will  be  excessive. 

If  stabilization  were  used,  the  permissible  errors  in  the  stabilized 
trajectory  would  be  the  same  as  those  given  in  the  preceding  section 
for  the  case  of  no  stabilization.  The  total  change  in  the  angular  orienta¬ 
tion  of  the  stabilized  platform  would  be  no  greater  than  the  scalar  sum 
of  the  angular  changes  about  three  mutually  perpendicular  axes.  If 
the  angular  rate,  Ap,  _a.bout  any  axis  did  not  exceed  one  third  of  the 
permissible  value  of  co  for  the  case  of  no  stabilization,  the  trajectory 
errors  due  to  these  rates  would  not  be  excessive: 

I  '^^permissible  ^  —  (V3)  I  “ without  stabilization^ 

where  Ap  includes  the  designed  motion  of  the  system  with  respect  to 
inertial  space  (if  any)  as  well  as  the  drift  rate.  Since  the  drift  rates 
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permitted  by  this  criterion  are  likely  to  be  well  within  the  capabilities 
of  available  instrumentation,  the  selection  of  stabilization  equipment 
might  be  dictated  by  considerations  other  than  accuracy  requirements, 
such  as  cost,  weight,  size,  etc. 

In  many  instances  a  stabilized  platform  can  be  used,  but  it  will 
not  be  possible  to  mount  all  the  instruments  on  the  platform.  In  this 
case  it  will  be  necessary  to  measure  and  record  the  orientation  of  the 
platform  with  respect  to  a  body-fixed  coordinate  system,  and  it  may 
also  be  necessary  to  measure  the  bending  angles  for  each  data  station. 

If  relative  velocities  or  accelerations  are  to  be  measured  by  non- 
inertial  instruments,  rate  gyros  or  similar  instruments  could  be  used 
to  measure  the  angular  velocity  components  (at  each  data  station  unless 
it  can  be  shown  that  binding  will  not  contribute  significantly  to  these 
rates).  Alternately,  R  can  be  determined  with  respect  to  the  body 
axes  by  numerical  integration,  and  it  will  ni^  be  neces_sary  to  measure 
oi  .  The  latter  method  is  preferable  unless  and/or  a^  are  wanted. 

In  view  of  the  criteria  given  in  the  last  section,  it  appears  that  the 
angles  should  be  recorded  to  the  nearest  milliradian;  the  other  accuracy 
requirements  are  the  same  as  those  considered  in  the  last  paragraph. 

It  should  be  mentioned  that  some  vehicles  will  use  a  stabilized 
system  which  has  one  axis  "vertical"  and  one  axis  pointing  either  in 
the  direction  of  travel  or  to  geodetic  north.  The  word  vertical  is 
used  loosely  here;  an  exact  definition  cannot  be  given  without  consider¬ 
ing  a  specific  application  and  a  specific  method  of  determining  the 
vertical.  A  stabilized  system  of  this  type  is  useful  in  certain  types  of 
satellites  and  missiles  and  for  navigation  purposes  in  aircraft.  The 
system  will  rotate  slowly  as  the  vehicle  changes  its  position  above  the 
earth,  but  the  rotation  is  negligible  for  present  purposes. 

An  inertial  (non-rotating)  system  is  useful  for  deep-space 
vehicles,  for  certain  types  of  satellites,  and  for  ICBM's.  The 
orientation  of  such  a  system  with  respect  to  the  earth  at  the  beginning 
of  the  scoring  encounter  can  be  computed  with  sufficient  accuracy  given 
the  initial  orientation,  the  initial  location,  and  the  location  of  the 
intercept. 

The  last  method  of  correcting  for  angular  motions  is  to  use 
instruments  to  measure  these  motions  (including  bending)  and  then  to 
take  them  into  account  in  data  reduction.  A  minimum  of  three  angular 
rate  gyros,  sensing  angular  rates  about  each  of  three  mutually  perpen¬ 
dicular  axes  in  the  instrumented  object,  would  be  sufficient  to  establish 
the  changes  in  orientation  which  occur  during  the  course  of  an  attack. 

The  permissible  errors  for  these  rate  gyros  would  be  the  same  as 
those  for  the  gyros  used  for  stabilization,  as  given  above.  The  accuracy 
requirements  on  this  instrumentation  are  thus  not  critical.  Rate  gyros 
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are  commonly  used  for  measurements  of  this  type,  but  other  instru¬ 
ments  are  available.  The  angular  orientation  could  also  be  obtained 
from  photographs  of  the  object,  but  it  is  believed  that  this  procedure 
is  not  generally  applicable  and  might  prove  troublesome  in  data 
reduction. 

Rate  gyros  have  long  been  used  to  determine  orientation  as  a 
function  of  time  in  tests  made  on  rockets  and  aircraft.  Reports  giving 
the  details  of  this  use  are  readily  available. 

Integrating  the  angular  velocity  components  to  obtain  the  changes 
in  the  orientation  angles  is  a  troublesome  process  if  the  time  from 
initial  conditions  to  the  end  of  the  scoring  encounter  is  at  all  long.  The 
problem  of  obtaining  initial  conditions  will  usually  limit  the  application 
of  this  method  to  vehicles  such  as  aircraft  which  characteristically 
have  small  motions  about  a  nominal  orientation  v/hich  can  normally  be 
used  as  the  initial  condition  without  introducing  excessive  error.  It 
can  be  shown  that  the  angles  obtained  by  simple  integration  of  the 
angular  rates  do  not  have  an  inherent  rotational  order  and  as  a  result 
do  not  define  unique  orientations  of  the  vehicle;  the  discrepancies 
will  be  small,  however,  if  the  angles  are  small.  The  discrepancies 
can  theoretically  be  eliminated  by  using  a  simultaneous  integration 
process  employing  appropriate  coordinate  transformations,  but  the 
supposed  improvement  may  be  lost  as  a  result  of  roughness  in  the 
data  or  of  errors  stemming  from  the  vastly  increased  complexity  of 
the  calculations.  In  general,  the  probability  of  obtaining  good  accuracy 
over  a  long  period  of  integration  is  low.  In  short,  the  method  can 
hardly  be  recommended  except  as  a  last  resort. 
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APPENDIX  3 

RATE  CHARACTERISTICS  OF  LINEAR  RELATIVE  TRAJECTORIES 


Introduction 


In  problems  of  scoring  missile-target  encounters,  the  angular 
rates  of  one  of  the  objects  as  viewed  from  the  other  may  be  important 
variables.  If  the  measuring  equipment  is  sensitive  to  the  angular 
rates,  for  example,  the  errors  of  observation  will  be  dependent  upon 
these  angular  rates.  For  planning,  it  is  desirable  to  know  at  least 
the  angular  rates  and  the  ranges  of  variation  of  the  angular  coordinates 
in  order  to  be  certain  that  the  scoring  requirements  are  compatible 
with  the  allowable  scanning  rates  and  fields  of  view  of  the  sensors. 

It  may  also  be  useful  to  have  some  information  on  the  general  behavior 
of  the  angles  and  angular  rates  as  functions  of  the  time.  In  this  paper 
a  brief  investigation  of  such  matters  is  presented  without  regard  to 
the  particular  kind  of  equipment  to  be  used  for  scoring. 


Azimuth  and  Elevation  Angular  Rates 

The  angular  velocity  co  of  the  relative  range  vector  R  is  given  by 


R  X  R 


=  i(0  + 

X 


iio  + 
J  y 


kco 


(1) 


where  i,  j,  k  are  umt  vectors  in  an  X,  Y^Z  rectangular  Cartesian 
coordinate  system;  co  is  normal  to  both  R  and  R  . 


For  linear  relative  trajectoris,  R  is  given  by 

R  =  R  +  R  t 
o  o 


(2) 


Most  relative  trajectories  of  interest  are  well  represented  by  this 
linear  approximation.  Substitution  of  equation  (2)  into  equation  (1) 
gives 


CO 


R  X  R 

— =  M/R^ 


R 


(3) 


the  vector  M  being  a  constant  vector .  Clearly  |co  J  has  its  maximum 

when  I  R  I  is  a  minimum;  at  this  point  R  and  R  are  orthogonal  to 
each  other,  so  that 


'Max 


Relative  velocity 
Miss  distance 


(4) 


A-42 


Finite  values  of  |  cj  |  then  result  for  all  except  direct  hits. 

Certain  angular  rates  of  R  are  given  by  the  vectors  A  and  E, 
or  azimuth  and  elevation  angular  rates,  respectively  (see  Figure  1). 


Figure  1.  Azimuth  and  Elevation  Angles 

The  vector  E  is  given  by  ^ 

E  =  E  (i  cos  A  -  j  sin  A)  =  (co^  cos  A  -  sin  A)  E,^  "  (5) 

where  E ,  is  a  unit  vector  in  the  E  direction.  The  vector  A  is  given 
by 

*  • 

A  =  -kA  (6) 

which  is  not  necessarily  normal  to  R  .  The  vector  A  will  usually 
have  a  component,  positive  or  negative,  in  the  direction  of  R. 

• 

It  is  desired  to  be  able  to  put  bounds  on  the  angular  rates  A 

4  » 

and  E,  possibly  by  the  use  of  w.  Since  A  and  E  are  normal  to  each 
other,  and  the  magnitude  of  their  vectojr  sum  is  greater  than  or  equal 
to  the  length  of  co  ,  which  is  normal  to  R, 

co^  sA^  +  E"^  (7) 

This  puts  a  lower  bound  on  the  variables.  Also  since  E  is  normal  to 
R,  and  its  length  is  equal  to  the  projection  of  oo  in  the  direction  of  E  , 


(8) 
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As  I  (0  I  is  bounded  above  by  equation  (4),  this  same  bound  serves  as 

.  *  • 

an  upper  bound  for  the  magnitude  of  E.  Unfortunately,  |  A  |  is  not 
similarly  bounded. 

Putting  bounds  on  A  is  somewhat  more  troublesome.  One  has 


tan  A  =  X/Y 

A  =  (YX  -  XY)/(X^  +  Y^)  (9) 

For  linear  relative  trajectories,  this  reduces  to 

A  =  K/(X^  +  Y^)  (10) 


2  2 

which  has  its  maximum  magnitude  when  (X  +  Y  )  is  a  minimum.  For 
short-duration  trajectories,  the  maximum  magnitude  of  the  azimuth 
angular  rate  then  occurs  near  the  point  where  the  relative  coordinates 
2  2 

satisfy  (X  +  Y  )  =  minimum.  Ordinarily,  this  point  is  different  from 
the  minimum-miss-distance  point.  Clearly 

|A|j^ax-  |K  |/miss  distance^  (11) 


To  go  further  with  bounding  |  A  | ,  an  expression  for  A  in  terms 
of  w  will  be  derived.  One  has  the  fact  that  the  sum  of  the  projections 

of  A  and  E  on  oj  is  equal  to  the  length  of  w  ,  or 


«  • 

(A  +  E  )  •  Z 


2 

0) 


(12) 


Making  use  of  equations  (5)  and  (6)  produces 


'  *  2  ? 

-Ao)  +  E  =  o) 
z 

*  2*2 

-Aw  =  w  -  E 
z 


2  '2 

Equation  (13)  is  exact.  Since  w  >:  E  ,  one  has 


(13) 


(14) 
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From  eq^uation  (1)  one  has  for  linear  relative  trajectories 
u  =  (X  Y  -  Y  X  )/R^  =  K/R^ 

2  OO-  OO  ' 


(15) 


For  linear  relative  trajectories  one  may  substitute  from  equation 
(15)  into  equation  (14)  to  obtain 


A  < 


2 

R  CO 


X  Y  -  Y  X 
'  o  o  o  o ' 


(Relative  velocity) 
K 


(16) 


the  constant  K  being  given  by  K  -  I  ~  ^ 


Also,  since  j  E  |  is  the  projection  of  co  along  E.^,  then  |  E  |  is 
less  than  the  length  of  the  projection  of  w  onto  the  XY  plane,  or 


*  2  2  2 
E  <  w  +  w 
X  y 


(17) 


Substitution  of  equation  (17)  into  equation  (13)  leads  to 

I  A  I  >  I  CO  J  =  K/r2 


(18) 


Another  relation  comparable  to  equatio:^(  1  3)  may  be  derived  by 
considering  the  differential  displacements  of  R  due  to  the  angular 

motions.  In  Figure  2,  a  vector  R  undergoes  small  angular  displace- 

•  • 

ments,  A  At  and  E  At,  causing  a  rotation  of  amount  coAt. 

A  vector  R,  having  originally  the  direction  of  OT  in  the  X,  Y,  Z 

coordinate  system  of  Figure  2,  undergoes  small  angular  displace- 

•  • 

ments  due  to  A  and  E  acting  over  a  differential  time  At,  an_d  is  dis¬ 
placed  to  the  direction  of  OQ.  Consider  the  projections  of  R  onto  a 
sphere  of  radius  S  centered  at  O.  Triangle  TPQ  is  a  plane  right 
triangle  up  to  first-order  effects,  so  that 

(TP)^  +  (W)^  =  (TQ)^ 


(-A  cos  E)^  +  E^  (S  At)^  =  J-  (S  At)^ 


(19) 
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Hence 

cos'^E  +  (20) 

On  comparing  equation  (20)  with  equation  (13),  one  sees  that 

oj  =  -A  cos^E  (21) 

z 


Equations  (20)  and  (21)  are  applicable  in  the  general  case,  They 
may  also  be  derived  by  differentiation  (use  being  made  of  relations 
for  the  relative  rectangular  coordinates  in  terms  of  the  relative  polar 
coordinates  and  equation  (I))-.  Incidental  to  this  process,  one  may 

•  • 

obtain  relations  for  co  and  co  in  terms  of  A,  E,  A,  and  E. 

X  y  >  ,  , 

In  view  of  equation  (12)  an  immediate  consequence  of  equation  (20) 
is 


cos  E 


A 


cos  E 


(22) 


An  Example  of  a  Linear  Relative  Trajectory 


An  example  of  an  attack  with  linear  relative  velocities  is  the 
following: 
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at  t  =  0,  R  =  (50j  +  2000k)  ft 

* 

=  (40i  -  160j  -  4000k)ft/sec 

For  this  example  the  relative  rectangular  coordinates,  X,Y,Z,  are 
given  by 

X  =  40t 
Y  =  SO  -  160t 
Z  =  2000  -  4000t 

The  azimuth  and  elevation  angles  are  given  by 
tan  A  =  40t/(50  -  160t) 

tan  E  =  (2000  -  4000t)/(27,200t^  -  16,  OOOt  +  2500)2 


The  angular  rate  oj  of  R  is  given  by 

-  120,  OOOi  +  80,  OOOj  -  2000k 

oj  =  - - — _  ^  ^ 

4.  002,  500  -  16,  016,  OOOt  +  16,  027, 200t‘' 
The  angular  rate  A  =  (YX  -  XY)/(X^  +  Y^) 

A  =  2000/(2500  -  16, OOOt  +  27,  200t^) 

B 

The  angular  rate  E  is  given  by 


E  = 


E 


Z  «/x^  +  Y^  -  Z 

Jx2  +  Y^' 


/R 


r 

(2000  -4000t) 


-4000  *y  2500  -  16,  OOOt  +  27,  200t 

(-8000+  27200t) 


>Jz500  -  16,  OOOt  +  27,200t^ 


/(4, 002, 500 


-  16, 016, OOOt  +  16, 027,200t^) 
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Also 


=  4, 002, 500  -  16, 016, OOOt  +  16, 027,2001^ 


For  a  linear  relative  trajectory,  the  minimum  range  occurs 

at  t  =  -(R  »  R  )/(R  •  R  ),  with 

'  o  o"  '  o  o' 


R 


2 

Min 


These  formulae  may  be  obtained  by  differentiation  of 

R  =  R  +  2R  .  R  t  +  R^t^ 
o  o  o  o 

In  this  example,  tj^  «  0.49965  and  for  this  value  of  t, 

^min 

R  .  =  36.  084  feet, 

min 

Figure  3  is  a  plot  of  the  relative  rectangular  coordinates  X,Y,Z 
for  this  example.  These  coordinates  are  of  course  linear  functions  of 
time.  Figure  4  is  a  plot  of  R  vs  t,  which  is  linear  except  injthe  close 
neighborhood  of  the  minimum  range  point.  In  this  example  R^  and 

R  are  fairly  closely  aligned. 

Figures  5  and  6  are  plots  of  the  azimuth  angle  A  vs  t  and  the, 
elevation  angle  E  vs  t  for  this  example.  The  angular  rates  A  and  E 
are  shown  in  Figures  7  and  8.  Note  that  for  this  example  the 
maximum  angular  rate  A  does  not  occur  near  the  point  of  closest 
approach,  while  E  has  its  maximum  at  the  point  of  closest  approach. 

The  Angular  Rates  of  R  for  a  Certain  System  of  Angular  Coordinates 


One  system  of  angular  coordinates  measures  the  angles  of  the 
projection  of  R  on  the  ZY  and  XZ  planes  (the  angles  a  and  [3  of 
Figure  9).  The  angles  a  and  p  are  the  ones  measured  if  photo¬ 
potentiometers  are  used -as  sensors.  They  are  related  to  the  Cartesian 
X,Y,Z  coordinates  according  to 

R^  =  X^  +  Y^  +  Z^ 
tan  a  =  Y/Z 
tan  p  =  X/Z 


(23) 
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Example  of  Azimuth  Angle,  A  vs  t 
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Figure  9.  Angular  Coordinates  of  R 


This  coordinate  transformation  has  a  double-valued  inverse, 
but  in  application  the  values  of  Z  are  always  taken  positive,  so  that 
the  inverse  transform  is  given  by 

7  2  - 

Z  =  R/{1  +  tan  a  +  tan  p)^ 

Y  =  R  tan  a/(l  +  tan^a  +  tan^p)^  (24) 

X  =  R  tan  p/(l  +  tan^a  +  tan^p)^ 


Differentiation  gives 

a  =  -(YZ  -  ZY)/{Y^  +  Z^) 
P  =  -(XZ  -  ZX)/(X^  +  Z^) 


(25) 


It  is  desired  to  put  bounds  on  these  rates  in  terms  of  appropriate 
variables  of  an  attack.  One  has 
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-  ^  R  X  R  =  J_  i(YZ  -  ZY)  +  j(ZX  -  XZ)  +  k(XY  -  YX) 

R‘^  L 

Clearly,  one  has 

j(YZ  -  ZY)|  <  |w  I  R^ 

){XZ  -  ZX)  j  s  jo)  I  R^ 

Also,  one  has 

joj  I  <  Relative  velocity/R  :<  Relative  velocity /Miss  distance  (27) 

Then 


(26) 


=  |(YZ  -  ZY)/(Y^  +  Z^)|  < 


w  R 


|xz  -  ZX \/{yr  +  z^)  < 


Y  +  Z 

_MrL 

x^  +  z^ 


(28) 


For  attacks  with  constant  relative  velocities  all  the  quantities 
jYZ  -  ZY|,  |XZ  -  ZXj,  and  |w|R^ 

are  constants.  In  particular,  evaluating  w  at  the  minimum  range  point. 


“  ' ^ U 1 di^ tan^e^  ’  ~  (Relative  Velocity)(Miss  distance)^ 


and 


j«j  (Relative  velocity)(Miss  distance) 


*  I  ( Relative  velocity )(Miss  distance) 

P  I  -  - 2 - 1 - 

X  +  Z 


(30) 
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♦  2  2  % 

Since  |a  |  =  K./(Y  +  Z  ),  then  |a  J  has  its  maximum  if 

2  2^  *  22 
(Y  +  Z  )  is  minimum.  Similarly  |p  |  has  its  maximum  if  (X  +  Z  ) 

is  minimum.  These  points  in  general  are  different  from  each  other, 

and  are  in  general  different  from  the  minimum  range  point. 

Also,  (a  —  K.j/(Miss  distance)^,  with  a  similar  inequality 

for  Ip  Closer  bounds  are  given  by 

|i|  >  K^/R^ 

Ip  I  ^ 

The  components  oo  ,  co  ,  and  to  of  to  may  be  expressed  in  terms 

,  »  ^  y  ^ 

of  a,  p,  a,  and  p  by  means  of  the  coordinate  relations  (and  their 
derivatives)  given  by  equations  (24)  and  (25),  together  with  equation  (1) 
for  (0.  However,  no  important  or  obviously  useful  information 
appears  to  be  revealed  by  use  of  such  a  process.  The  bounds  that 
have  been  given  for  the  angular  rates  a  and  ^  appear  to  be  the  most 
useful  ones  for  planning. 
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APPENDIX  4 

MINIMUM  DATA  REQUIREMENTS  AND  VARIANCE  CONSIDERATIONS 
FOR  ONE-,  TWO-,  AND  THREE-STATION  MEASUREMENTS 

Introduction 


The  determination  of  the  position  of  an  inaccessible  point  by 
triangulation  is  one  of  the  oldest  methods  of  indirect  measurement. 

This  method  has  received  considerable  theoretical  attention  in  recent 
years,  particularly  in  error  analysis  studies.  Most  of  these  theoretical 
treatments  have  considered  angular  measurements  only  and  have  been 
concerned  with  the  over-determination  of  the  point  so  that  the  errors  in 
the  indirectly  measured  quantities  can  be  reduced  or  minimized  by  use 
of  weighting  factors  or  by  making  a  least- squares  adjustment.  The 
study  reported  here  was  not  concerned  with  the  minimum  sets  of  data 
required  for  one-,  two-,  and  three- station  measurement  systems 
and  was  not  restricted  to  angular  measurements.  The  usual  angle-only 
triangulation  methods  were  of  course  included. 

In  rectangular  coordinates,  the  minimum  information  is  a  set 
of  values  for  the  (X,Y,Z)  coordinates  of  the  point  in  question.  In 
polar  coordinates,  the  data  required  would  be  the  distance  from  the 
origin  to  the  point  plus  the  azimuth  and  elevation  angles  of  the  line  of 
sight  (R,A,E).  If  measurements  are  taken  from  more  than  one 
position  in  the  coordinate  system,  not  all  of  the  quantities  (X,Y,Z),  or 
equivalently  (R,  A,E),  need  be  measured  from  each  position..  For 
example,  an  azimuth  angle  measured  at  one  station,  an  azimuth  and  an 
elevation  angle  measured  at  a  second  station,  and  the  distance  between 
the  two  stations  are  sufficient  to  establish  the  position  of  a  point  P 
if  the  point  does  not  fall  on  the  extended  line  segment  which  connects 
the  two  measuring  stations. 

A  few  general  remarks  are  made  concerning  one-station 
measurements,  while  all  possible  combinations  of  minimum  sets  of 
data  are  exhausted  in  a  systematic  manner  for  two-station  measure¬ 
ments.  Only  a  fev/  selected  cases  of  minimum  data  are  considered 
for  three-station  measurements.  The  results  are  presented  in  an 
abbreviated  form  for  easy  reference.  Equations  and  procedures  are 
developed  for  specific  cases  to  show  how  the  variances  (mean  square 
errors)  in  the  computed  results  are  related  to  the  variances  in  the 
measured  data  and  how  either  the  ranges  of  variation  or  the  variances 
of  the  measured  quantities  are  limited  by  the  requirement  that  the 
variance  in  a  given  result  should  not  exceed  a  specified  value.  A 
method  is  given  for  weighting  the  measured  quantities  so  that  the 
variance  of  a  computed  quantity  can  be  minimized  if  redundant  data 
are  available. 
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One-Station  Measurements 


Figure  1  shows  a  point  P  and  a  single  measuring  station  at  point 
O  which  is  taken  as  the  origin  of  the  X,  Y,  Z  coordinate  system.  The 
position  of  point  P  can  be  determined  by  measuring  its  Cartesian 
coordinates  (X,y,Z),  but  it  will  usually  be  more  convenient  (or  even 
necessary)  to  measure  the  spherical-polar  coordinates  (B,,A,E).  If 
either  of  these  sets  is  specified,  the  other  set  is  easily  obtained 
through  the  simple  relations 


Z 


Figure  1 


One-Station  Measurements 


X  =  R  cos  E  cos  A 
Y  =  R  cos  E  sin  A 
Z  =  R  sin  E 


2  2  2  - 
R  =  (X  +  Y  +  Z  )2 

A  =  tan"-^(Y/X) 

E  =  sin'^(Z/R) 


If  the  differentials  are  regarded  as  small  errors  or  deviations 
from  a  correct  set,  the  approximate  errors  in  rectangular  Cartesian 
coordinates  and  in  spherical-polar  coordinates  are  represented  by  the 
differentials  as  follows: 

dX  =  dR  cos  E  cos  A  +  R  [  -sin  E  cos  AdE  -  cos  E  sin  AdA] 
dY  =  dR  cos  E  sin  A  +  R  [  -sin  E  sin  AdE  +  cos  E  cos  AdA] 
dZ  =  dR  sin  E  +  R  cos  EdE 


and 


A-60 


dR  = 


(X^  +  +  Z^)2 


[XdX  +  YdY  +  ZdZ] 


dA  = 


2  2 
(X  +  Y  ) 


XdY  -  YdX] 


dE 


(X--  +  Y‘') 


2.4  L 


dZ  - 


Z(XdX  +  YdY  +  ZdZ) 

?  2  ?  - 
(X  +  Y  +  Z  )* 


No  redundancy  of  data  can  result  when  only  one -station  measurements 
are  used. 


Two-Station  Measurements 
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The  point  O,  located  midway  between  the  two  measuring 
stations,  and  M2,  will  be  the  point  to  which  the  position  of  the 

observed  point,  P,  is  to  be  related.  The  azimuth  and  elevation  angles 
at  point  O  will  be  designated  as  a  and  e  ,  respectively. 

The  separation  of  the  two  stations  will  be  designated  as  f  .  It 
will  be  assumed  that  this  distance  has  been  determined  exactly,  previous 
to  the  measurements  made  on  the  point  P.  Though  this  is  not  strictly 
true,  the  error  in  measuring  I  should  certainly  be  considerably  less 
than  the  error  in  measuring  or  p2.  Following  this  same  line  of 

reasoning,  it  is  assumed  that  point  O  is  located  at  the  exact  midpoint 
of  f  . 


To  determine  all  possible  minimum  sets  of  data  from  the  two 
measuring  stations,  we  observe  that  the  number  of  combinations  which 


can  be  made  with  six  quantities  when  taken  three  at  a  time  is 


6; 

■3TC&"-  T)  '■ 


or  20  different  minimum  sets.  Of  the  20  possible  sets,  some  may  be 
discarded  since  they  do  not  afford  sufficient  information.  As  an  ex¬ 
ample,  if  E^,  pj^,  and  A2  are  the  three  measured  quantities  being  con¬ 


sidered,  the  point  P  is  not  uniquely  determined  since  there  will  be  two 
points  on  the  elevation  cone  of  station  Mj^  which  are  located  a  distance 


p^  from  M^  and  fall  in  the  plane  established  by  A^.  This  condition  is 
pointed  out  in  Figure  3.  Both  points  P^  and  P2  satisfy  the  same 
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Because  of  the  symmetry  between  the  two  measuring  stations, 
some  of  the  cases  may  be  obtained  as  obvious  extensions  of  others.  An 
example  of  this  symmetry  effect  is  afforded  by  the  two  minimum  data 
sets  (Aj^,  E^,  A^)  and  (A^,  E^,  Aj^).  The  equations  for  the  first  set 

can  be  transformed  to  the  equations  for  the  second  set,  and  vice  versa, 
by  interchanging  subscripts  i  and  2. 

Still  other  cases  which  may  be  discarded  are  those  which  turn 
out  to  be  one-station  sets  of  measured  quantities.  Only  two  such  cases 
exist,  namely  (Aj^,  Ej^,  pj^)  and  (A^,  E^,  P2)- 

The  cases  which  will  not  be  considered  in  detail,  either  because 
the  point  P  is  not  uniquely  determined  or  because  they  are  one-station 
measurements,  are  as  follows; 


1. 

(A^, 

^2- 

Pi) 

2. 

(A^, 

E^, 

P2^ 

3. 

El, 

P2> 

4. 

(A^. 

El, 

Pi) 

5. 

(El, 

^2* 

Pi) 

6. 

(Ep 

E^, 

P2> 

7. 

(E^. 

Pi' 

P2> 

8. 

(^2- 

Pr 

P2^ 

9. 

(A^, 

El, 

Pi) 

10. 

(A^, 

^2' 

P2> 

The  equations  for  five  out  of  the  remaining  ten  cases  will  be  given  in 
the  following  paragraphs;  the  results  for  the  other  five  cases  are 
readily  obtained  from  symmetry  considerations .  In  these- equations,  angles 
a,  Aj^  and  A^  all  have  the  same  sign;  the  quadrant  of  a  will  have  to  be 

determined  from  the  inequality  Ia^I  £  |a|  |A^  |.  Note  that  only  two 

out  of  the  five  cases  are  completely  determinate;  however,  the 
indeterminate  cases  can  sometimes  be  solved  by  restricting  the  fields 
of  viev/  of  the  measuring  devices. 
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Case  I 


where 


Case 


.  Measured  Quantities:  i,  A^,  E^,  A^, 
G  +  H  -  J  •  2 


S.  =  i 


G  +  H  -  J  2 
1 - ^ - 1 


.  -1 
=  sin  , 


.  -1 
a  =  sin 


[  G  +  H  -  4J] 


MN 


[  G  +  HN  -  4J] 


2  2 

G  =  cos  sin  (A2  -  Aj^) 


H  =  4  sin  A- 


J  =  cos  AjSin  A^  cos'^E^sin  (A^  -  A^) 

K  =  4G 

L  =  2  sin  A^  sin  E^ 

M  =  2  sin  A^  sin 
N  =  cos  Ej^ 

II.  Measured  Quantities:  i  ,  Aj^,  E^,  E^  (valid  when  E  >  E^) 

n  -  G^JK 

^  ,2  4H 


€  =  Sin 


G^JK  ,  L  \i 

“J 

L  [4JG‘^K  +  HL] 


r 


a  =  sin  tan  • 


2G^K 


|_  '  [4JG^  +  HL  -  M^G^K^] 
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where 


G  = 
H  = 
J  = 
K  = 

L  = 

M  = 


cos  Aj^  cos  Ej^  sin  E^ 

sin^E^  - 

sin  Ej^ 

sin^E^  + 

sin^E 

H 

o 

NJ. 

•3  ■  2„ 
-3  sxn  E 

2 

^  -  sin  Ej 

2  sin  Ej 

sinE^ 

Case  III.  Measured  Quantities:  i,  p^, 


R 


=  { 


+  i  +  iHJ 


<  =  cos  < 


a  =  cos  < 


[  J^cos^E^  +  +  iHJL] 

[p^  +  i^G  +  iHJ]  ^ 


[p^H^  +  i  -  H'^)}  +  iHJ{2H'^-  1)] 


,2f  1 


[p^  M  +  -  M  +  H^L}  +  iHJL] 


where  e  is  positive  if  is  positive,  and 
G  =  (-  -I"  +  cos^A^  cos^  E^) 

H  =  cos  A^ cos  E^ 

J  =  [p^  -  i^(l  -  cos^A^.cos^E^)]  ^ 

12  2  2 
K  =  (  -  sin  E^cos  E^cos  A^) 

2  2 
L  =  cos  E^  -  sin  E^ 

M  =  cos^E^ 
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Case  IV.  Measured  Quantities;  i,  pj^,  p^, 


R  =  ^  (2G  - 


-R 

e  =  cos  ^ 


i  [H^  -h  J^(H  +  i^)^]  ^ 
[2G  -  H] 


a  =  tan 


where  the  sign  of  «  is  indeterminate,  and 

^  2^2 

G  =  +  P2 

„  2  2 
H=  Pi  -  P2 


J  =  tan  A, 


Case  V.  Measured  Quantities:  f  ,  p^,  Aj, 
G  -  4H' 


R 


_  r  2  ,  ,  2  r  G  -  4Hl  1  ^ 


•  -1  , 
«  •-  sm  J 


a  =  sin 


'{ 


7  7  7  7  -  ”> 

2[pJg^ 


[4piG^  +  f  ^G(G  -  4H)] 


2K 


[G^  -  4{GH  -  J^)]  ’ 


i^} 


where  the  sign  of  €  is  indeterminate,  and 
G  =  sin  (A^  -  A^) 

H  =  cos  AjSin  A^ 

J  =  sin  A^ 

K  =  sin  A^  sin  A^ 
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Expressions  for  the  differential  errors  in  the  computed  quantities 
R,  a,  and  c  resulting  from  small  (differential)  errors  in  the 
measured  quantities  p^,  p^.  A^,  A^,  E^,  and  E^  have  been  obtained 

by  the  straightforward  method  of  taking  differentials.  While  the  method 
is  relatively  simple,  the  resulting  expressions  are  exceedingly 
cumbersome  and  have  not  been  included  here.  For  possible  future 
computational  work  and  for  reference,  these  differential  expressions 
will  be  kept  on  file. 

These  differential  error  formulas  may  be  used  to  set  limits  on 
the  allowable  errors  of  the  measuring  devices  if  the  accuracy  required 
for  the  computed  quantities  is  specified.  To  facilitate  this  discussion 
of  allowable  errors,  two-station  geometry  will  be  considered  in  which 
the  azimuth  angle  is  measured  at  one  station  and  both  the  azimuth  and 
elevation  angles  are  measured  at  the  other  station;  this  is  one  of  the 
two  completely  determinate  cases  considered  above.  For  simplicity, 
the  computed  distance,  R,  is  related  to  the  station  at  which  both 
azimuth  and  elevation  angles  are  measured.  The  expression  for  R  in 
terms  of  the  measured  angles  can  be  written  as  follows: 


JR 

I 


sin  A- 

- ^ -  =  f(A.,A^,E,) 

cos  Ej  sin{A2  -  A.^) 


where  f  is  the  station  separation.  The  variance  (or  average  square  of 
the  error)  in  R  is  defined  by 


Limit 
Ni — f 


2 

where  R  =  R(X  ,  X^,  ....  Xj^)  and  is  the  variance  in  variable  X^. 

n 

The  result  obtained  for  this  case  can  be  written  as  follows; 


where 


sin  A^  cos(A2  " 

_ 

cos  Ej^  sin  (A^-Ap 


2 


A-67 


sin(A2  -  A^)  cos  A^  -  sin  A^  cos  {A^  -  A^) 

- ~~1 - 

cos  E  j^  sin  (A^  " 

sin  A^  sin  Ej^  ^ 

cos^  E|  sin  (A^  -  A^^) 


For  convenience,  it  is  now  assumed  that  the  variances  in  the  angular 
measurements  are  all  equal;  i.e.  , 


2 

‘^A 


1 


Z 

m  0- 


and  the  expression  for  the  variance  in  R  may  then  be  rearranged  to 
give 


1 _ 

f‘^(Aj,A2,Ep 


} 


The  right-hand  side  of  this  expression  is  dimensionless  and  plots  of 

this  expression  and  of  the  expression  for  (f/R)^  may  be  obtained  as 
functions  of  the  angles  Aj^,  Ej,  and  A^.  Such  plots  are  given  in 

Figures  4  through  15  which  follow.  The  parameter  A^  is  used  as  the 

connecting  variable  between  the  two  sets  of  curves. 


For  showing  how  these  curves  may  be  interpreted,  we  consider 
the  following  design  problem.  Suppose  the  variance  in  R  is  required 

to  be  less  than  tr?,  over  a  range  of  values  for  R,  say  R  <R.  , 

The  fields  of  view  of  the  angle-measuring  devices  and  the  variances 

(T^  (for  all  angles  measured  by  these  devices)  are  specified.  The 
question  to  be  asked,  then,  is  how  long  must  the  base  line  be  for  all 
these  requirements  on  our  detecting  system  to  be  fulfilled- -assuming 
that  they  can  be  fulfilled.  The  first  step  is  to  determine  values  for 

A^  by  using  the  curves  of  (  ^  j  vs  A^  for  all  E^  —^1 

“^Min  ^ 

2  2 

given  values  of  cTj^  ,  cr  (all  angles),  A^^ 


and  R 


Min 


Max' 


Max 

Some  of  the 


values  of  A_ 


'Min 


that  would  otherwise  be  acceptable  may  not  be  within 


(degrees ) 


(degrees) 


(degrees) 


(degrees) 
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(Elevation  Angle  E  .  =  ±45 


A-75 


(degrees) 


(Elevation  Angle  E.  =±60 
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ID 

rj 

(\l 

• 

* 

• 

O 

o 

d 

Pd 

O 

O 


(degrees) 


(degrees) 


o 

00 

sO 

fsj 

r\J 

o 

—4 

o 

o 

O 

O 

d 

d 

d 

Dc; 

d 

d 

A-80 


the  field  of  view  of  the  device;  in  such  cases  the  field  limit  should  be 

2 

used.  The  values  of  {f /R)  corresponding  to  each  set  (A.  ,  A^ 

2  Xfin  Min 

E.)  can  now  be  read  from  the  (f/R)  -vs-A-  curves.  Since  R.  ,  is 
1  '  '  '  2  Max 

known,  the  required  value  of  I  may  be  readily  computed  from  the 
2 

largest  value  of  (f/R)  determined  in  this  manner. 

Two-Station  Angle-Only  Measurements--One  Angle  at  Each  Station 


It  is  desired  to  investigate  the  geometrical  and  error  relations 
for  ranging  with  two-station  measurements  of  angles  only.  In  Figure  16, 
the  observed  object  is  at  C,  the  base  line  AB  is  of  known  length,  D, 
and  angles  A  and  B  are  determined  from  the  azimuth  and  elevation 
angles  measured  at  A  and  B  by  the  relations  cos  A  =  cos  A^  cos  A^, 

and  cos  B  =  cos  B  .  cos  B_  where  A  .  and  A_  are  the  azimuth  and 
A.  ii*  A  Xj 

elevation  angles  measured  at  A,  etc. 

C 


Figure  16 

Two-Station  Scoring  Geometry 


From  the  geometry, 

Rj^/D  =  sin  B/sin(A  +  B) 
R2/D  =  sin  A/sin(A  +  B) 


•  2.  ,  .  2„ 
sin  A  +  sin  B 

sin^  (A  +  B) 
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The  last  equation  was  used  to  compute  R/D  as  a  function  of 

(A  +B)  for  B  held  fixed  at  each  of  several  values:  15°,  30°,  60°,  90°. 
The  results  are  plotted  in  Figures  17  through  20.  It  was  noted  that 

larger  values  of  R/D  occur  for  (A  +  B)  near  180°,  a  condition  that 
would  exist  when  long  ranges  are  involved.  Additional  computations 
were  made  to  cover  the  cases  with  a  subtended  angle  C  on  the  interval 

0-5°,  and  Figures  21  through  24  show  plots  of  R/D  vs  B  for  the 
subtended  angle  C  held  fast.  These  plots  show  that  for  a  constant 
angle  C,  maximum  range  occurs  when  R  is  normal  to  the  base  line,  AB. 
Also,  R/D  is  a  highly  non-linear  function  of  C  as  C  approaches  zero. 

The  relations  connecting  the  angular  errors  made  in  measuring 
the  angles  A  and  B,  the  errors  of  R,  and  the  length,  D,  of  the  base 
line  AB  are  of  interest.  Plots  17  through  24  show  that  the  critical 
situations  for  accuracy  are  those  with  either  A  or  B  small  and  the 
other  large,  so  that  the  subtended  angle  C  is  small.  In  such  situations, 
R/D  is  large,  and 

R/D  ~Rj/D~R2/D 
Making  use  of 

R/D  ~  R^/D  =  sin  B/sin  (A  +  B), 

and  taking  differentials,  ignoring  for  the  present  any  errors  in  D,,  gives 

dR  cos  B  dB  sin  B  cos  (A  +  B)  (dA  +  dB) 

W  ~sm  (A  +  By  (a  +  B> 

Upon  setting  the  measurement  errors  in  A  and  B  equal  to  each  other, 
dA  =  dB ,  and  multiplying  dR/D  by  D/RdA,  there  results 


dR  _  cos  B 

RdA  "  R  sin  (A  +  B) 


2  sin  B  cos  (A  +  B) 
sin^(A  +  B  ) 


R 


H(A,B> 


Plots  of  H(A,B)  vs  (A  +  B)  for  B  held  fast  were  made  '(Figures 
25  through  27).  This  function  becomes  large  without  limit  for 

(A  +  B)  =  180°.  If  B  <90°,  and  (A  +  B)  £  90°,  the  two  terms  in  H(A,B) 
have  the  same  algebraic  sign  for  dA  =  dB.  As  (A  +  B)  tends  toward 

180°,  the  second  term  predominates,  sp  that 


R/D  vs 
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Figure  21 

R/D  vs  B  for  Subtended  Angle 


A-89 


o 

o 


Q 


vs  B  for 
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H(A,B)  vs  { 
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dR 

E3X 


R 


-2  sin  B  cos  (A  +  B) 


sin  (A  +  B) 


2/ sin  (A  +  B) 


If  C  S  6°  and  dA  «  C/2,  this  asymptotic  relation  gives  good  accuracy. 
One  has  also  in  this  neighborhood 


~  2/sin  B 


This  last  relation  may  be  useful  in  determining  the  relations 
among  the  values  for  angular  and  ranging  error,  the  range,  and  the 
length,  D,  of  the  base  line.  If  all  except  one  of  these  quantities  is 
known,  the  remaining  quantity  may  be  found.  Ratios  could  also  be 
found  for  more  than  one  unknown. 

As  an  example  of  the  use  of  the  last  equation,  suppose  it  is  given 

that  R  =  2,  000  ft,  B  =  15°,  dA  =  3  milliradians,  dR  =  100  ft.  What 
length,  D,  of  base  line  is  consistent  with  this  accuracy?  From  this 
equation  D  =  927.  3  ft  and  this  or  larger  values  of  D  would  give  the 
required  accuracy.  A  more  accurate  computation  shows  that  this  value 
of  D  corresponds  to  a  range  error  of  132.7  ft;  the  last  equation  is  thus 
a  fairly  good  approximation  even  for  angle  C  as  large  as  9  .It  should 

always  be  verified  that  2dA  «  180°  -  (A  +  B)  before  making  use  of  this 
last  relation.  As  another  example,  suppose  R  =  2000  ft,  D  =  100  ft, 

dR  =  100  ft,  B  =  15°.  What  size  errors  dA  are  permitted?  Using  the 
last  equation,  dA  =  0.  32  mils;  angular  errors  of  this  size  or  smaller 
would  give  range  errors  within  the  100-ft  limit. 

Problems  of  these  two  types  may  also  be  solved  by  use  of  the 
curves.  Figures  28  through  30  show  plots  of  R/D  and  dR/RdA  vs  (A+B) 
for  B  held  constant.  The  value  of  (A  +  B)  can  be  determined  from  one 
of  these  ratios,  R/D  or  dR/RdA,  and  the  other  ratio  may  be  read  from 
the  graph  at  the  corresponding  value  of  (A  +  B).  The  desired  quantity 
can  then  be  computed  using  the  value  of  the  ratio,  together  with  the 
values  of  the  remaining  variables. 

If  the  subtended  angle  C  is  not  small,  say  C  >  5°,  then  the 
approximation  R/D  ~  R^/D  is  no  longer  valid.  In  this  case  more 

accurate  expressions  should  be  used  for  the  differentials.  On  taking 
the  differential  of  the  R/D  equation  and  multiplying  by  D/RdA,  the 
result  for  dA  =  dB  is 
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Figure 
dR/RdA  vs  (A  +  B] 
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dA  _  1 
RdA  ~  2 


sin  B 

sin  (A  +  B) 


(-2  sin  B  cos  (A  +  B)  ^ 
sin^(A  +  B ) 


cos  B 
sin  (A  +  B) 


,  sin  A  ( -2  sin  A  cos  (A  +  B)  ^  cos  A  \ 
sin  (A  +  B)  sin2(A  +  B)  (A  +  B)j^ 

This  relation  for  can  be  used  with  the  equation  for  (R/D) 

to  solve  problems  of  the  type  of  the  two  examples.  One  of  the  relations 
can  be  used  to  find  (A  +  B),  given  B,  and  this  value  of  (A  +  B)  can  then 
be  used  to  find  the  unknown  ratio  which  is  to  be  solved  for  desired 
quantity  (D,  dA,  or  dR). 


In  conclusion,  it  may  be  stated  that  the  most  severe  requirements 
are  placed  upon  the  accuracy  of  angles  or  the  size  of  D  when  the  sub¬ 
tended  angle  C  is  small  or  either  A  or  B  is  small.  Under  such 
conditions,  the  use  of  this  method  of  ranging  and  position  finding  by 
angle-only  measurements  is  impractical.  If  range  as  well  as  angle 
information  is  available,  the  accuracy  requirements  would  be  greatly 
reduced;  if  redundant  data  are  available,  least  squares  or  a  similar 
method  could  be  applied  to  improve  the  estimate  of  the  position  of  the 
observed  object. 


Also,  measurements  of  angles  A  and  B  alone  could  give  good 
ranging  and  position  data  for  short  ranges,  even  with  a  short  base  line,, 
the  subtended  angle  C  being  large.  Such  data  could  be  extrapolated  to 
give  a  trajectory.  Unless  the  angular  measurements  are  made  with 
extreme  accuracy,  it  appears  probable  that  at  extreme  ranges  of 
interest  for  scoring  the  extrapolated  trajectory  would  be  more  reliable 
then  the  trajectory  given  by  triangulation  at  these  ranges. 


Three-Station  Measurements 

The  three-measuring-station  geometry  is  as  indicated  in 
Figure  31. 

The  X^,  Yj,  X^,  ’^2'  ^3’  ^3  ^^®  pfa-ne  determined  by 

the  three  measuring  stations.  The  distances  f  and  i  ^  and  the  azimuth 
angle  A^,  are  fixed  previous  to  any  measurements  of  point  P,  and  the 
errors  in  the  magnitudes  of  i  ^  i  and  A^,  are  assumed  to  be  negligible 
in  comparison  to  the  errors  introduced  through  the  measurements  of 
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point  P.  Only  a  few  special  cases  will  be  considered  for  the  three- 
station  situation,  since  any  multiple -station  measurements  may  be 
reduced  ultimately  to  a  two-station  measurement. 

The  position  of  point  P,  for  all  cases,  will  be  referred  to  a 
coordinate  system  located  at  the  second  measuring  station,  i.e.  ,  to 
the  X^i  Y  2^,  7^2  reference  frame. 

Case  I.  Measured  Quantities:  f  ,  i  -A-pi  Rji  R2>  R3 
--  Measured  Directly 

"^2  = 
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E., 

=  cos 

-1 

1 

K 

2 

where  the 

sign 

of  I 

L 

=  2- 

G  = 

sin 

Ap 

H  = 

cos 

J  = 

£  ( 

-  ^7 

1  \ 

.  3 

2 

K  = 

^  2 

(rJ 

1  + 


M 


ZR^  2 


M  =  H  + 


(4) 


Case  II.  Measured  Quantities;  Ap,  A^  A^, 


R2  =  +  G^)  +  G^J  (J  -  2HK) 


>? 


A^  --  Measured  directly 


E  2  =  tan 


-  2HJK| 


where 


G  =  tan  E, 


H  =  i  j  sin  Aj 
J  =  2 

K  =  cos  (A2  -  Ap) 


Case  III.  Measured  Quantities:  i  A„,  A.,  A.,  E, 

I  ^  V  i  J  ^ 

+  2HJL  !  ^ 


=  +  2HJlJ>^ 

if  H[M  +  J]  1_ 

1  L  HN  +  i  ^IJJ  ' 


A2  =  tan 


A-lOO 


--  Measured  directly 

where 


G 

= 

f  sec  E2 

H 

= 

sin  Aj^ 

I 

= 

cos  A, 

J 

_  Cl  n  1  A  — 

f  ^  sm 

K 

= 

sin  (Aj  -  Aj^) 

L 

= 

cos  (A^  -  A^) 

M 

sin  A3 

N 

= 

cos  A3 

Case  IV.  Measured  Quantities:  ly  ly  Ap,  A^,  A^,  E^ 


^2  = 


H  (1  +  G  )  +  G  J  (J  +  2HK) 


Measured  directly 


E^  =  tcun 


■‘{ll 


+  2HJK] 


where 


G  =  tan  E^ 

H  =  i2  (^3  - 

J  =  i  J  sin  (A^  -  A^) 
K  =  cos  A^ 

Weighting  Redundant  Data 


The  accuracy  with  which  the  position  of  a  point  in  space  may  be 
determined  with  respect  to  some  arbitrary  reference  frame  is 
completely  dependent  upon  the  accuracy  of  the  measuring  devices 
being  used.  A  quantity  which  is  not  measured  directly,  but  which  is 
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^2  = 


-1 


L 


1  + 


M 


where  the  sign  of  is  indeterminate,  and 
G  =  sin 


H  =  cos  At 


l(' 


J  = 


) 


K  =  f. 


ZR^  I  2 


M  =  H  + 


(4) 


Case  II.  Measured  Quantities: 


1’ 


2’  "^F’  '^1’  ^Z’  ^3 


^2  =  T 


H  (1  +  G  )  +  G  J  (J  -  2HK) 
Measured  directly 


■? 


=  tan 


-  2HJK| 


where 


G  =  tan  E, 


H  =  f  J  sin  Aj^ 

J  =  f  2  ^^^2  " 

K  =  cos  (A^  -  Ap) 


Case  III.  Measured  Quantities:  i  , ,  f.,,  A_,  A,,  A,,  E^ 

I  c.  ji  1.  5  a 


R2  =  +  2HJL 


A^  =  tan 


•  ij  H[M  +  J] 

i  L  hn  +  f  jij  J  r 
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computed  from  several  measured  quantities,  will  also  be  subject  to 
the  errors  introduced  through  the  measuring  devices.  Assume  that  X 
is  a  quantity  which  is  to  be  computed  from  several  measured  quantities, 
Z^,  . Z^,  and  that  a  functional  relationship  between  X  and  the 

measured  Z 's  exists.  This  relationship  could  be  written  as 


X  =  X(Z,.  Z,,  .  .  .  ,  Z  ) 

'1  2  n 

If  the  errors  in  the  measured  quantities  are  sufficiently  small,  the 
error  in  the  computed  quantity,  X,  is  given  approximately  by  the 
linear  relation 


dX  = 


ex 

0Z  , 


4. 


^7  4.  4- 


dZ 


where  dZ^  is  the  error  in  the  measured  quantity  Z^.  Further,  if  the 

error  distribution  in  the  measured  quantities  is  such  that  the  greatest 
majority  of  the  errors  are  sufficiently  small  for  the  Taylor  expansion 
to  be  approximately  valid,  and  if  the  errors  are  independent,  then 


_1_ 

N 


or  the  variance. 


2 

“■x 


is 


2 

(  f  cr^  + 

f  8X 

J  ■ 

CTx  = 

•  •  VFZTj  ^Z 
^  n'  n 

If  the  dZ's  have  a  mean  of  zero,  the  computed  quantity  dX  will 
have  a  mean  of  zero  also.  No  information  is  required  about  the 
distribution  of  errors  to  get  a  mean  of  zero;  however,  for  obtaining  the 
approximate  average  square  of  error,  additional  information  about  the 
distribution  is  required. 

It  is  to  be  stressed  that  once  the  functional  relationships 
between  the  measured  quantities  and  the  quantities  to  be  computed 
have  been  established,  the  differentials  of  these  relationships  may  be 
used  to  establish  limiting  conditions  on  the  variances  that  can  be 
tolerated.  If,  for  example,  in  the  above  relationship,  X  =  X(Z.j^,  Z^,-.., 

Z  ),  the  average  square  values  of  the  errors  involved  in  measuring  the 

Zj.  are  known,  then  the  approximate  variance  or  average  square  value 

of  error  in  the  computed  quantity  X  is  obtained  from 
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2 

(  f  2 

! 

ax  \2  2 

“■x  " 

\dZ^)  ^ 

}  '^2  ■ 

•  •  ^  V 

az  / 

n  ^ 

The  measurements  of  the  Z^'s  are  assumed  to  be  uncorrelated,  so  that 

for  N  measurements 

_  N 

(dZjdZ.)  =  ^  ^  (dZidZj)^  »  0 
k=l 
2 

The  coefficients  in  the  equation  will  in  general  turn  out  to  be  n 

rather  different  functions  of  the  dependent  variables.  If  the  variances 
2 

cr^  are  equal,  only  the  maximum  sum  of  the  coefficients  will  be  needed 

for  error-study  purposes.  In  general,  however,  it  will  be  necessary 
to  compute  several  "extreme"  sets  of  coefficients.  If  the  coefficients 
are  bounded,  they  can  be  used  to  set  limits  on  the  variances  such  that 

a  specified  will  not  be  exceeded;  unless  an  arbitrary  distribution 

of  the  allowable  error  is  made,  a  trial-and-error  process,  with  some 
attention  given  to  the  practical^limits  of  the  equipment,  will  be  needed 
in  general.  Of  course  if  the  <rf's  are  known  to  have  particular  values, 

then  the  maximum  {or  minimum)  variance  of  the  computed  quantity,  X, 
can  be  established. 


Some  comments  on  the  theorem  by  Tschebyscheff  might  be  of 
interest.  This  theorem  states  that  "the  probability  of  obtaining  a 
value  of  the  standardized  variable  which  is  numerically  less  than  or 


equal  to  a  specified  number  k  is  larger  than  1 


>  1 


"  That  is. 


Some  numerical  values  of  this  function  are  as  follows: 


1.  0 

1.41 

1.  58 

1.83  2.0 

2.  24 

3.  0 

3.  16 

10 

0 

.  50 

.  60 

.70  .75 

.  80 

.  889 

.90 

.  99 

k" 

From  this  table,  according  to  the  theorem,  it  can  be  seen  that 
independent  of  the  distribution  of  errors  and  independent  of  the 
magnitude  of  cr,  the  probability  of  measuring  X  to  within  say  1.  83  cr 

units  of  the  correct  value  X  is  greater  than  7  0  percent..  Thus  some 
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limits  on  the  reliability  of  the  measured  quantities  may  be  established. 

The  probability  of  an  error  not  being  greater  than  a  given  size 
may  be  stated  in  terms  of  the  variance,  according  to  Tschebyscheff 's 
theorem.  If  the  distribution  of  the  errors  is  known,  the  bounds  of 
these  probability  relations  may  be  made  more  precise.  For  the 
normal  distribution  the  probable  error  is  0.  6745cr.  This  statement  is 
to  be  interpreted  as  giving  the  size  of  the  error  such  that  one  half  the 
errors  are  less  than  that  size  and  the  other  half  are  greater.  According 
to  common  usage,  experimental  results  may  be  given  as  the  measured 
values  plus  or  minus  the  probable  errors. 

If  a  requirement  were  made  that  the  size  of  the  errors  be  less 
than  B  with  a  probability  C,  Tschebyscheff 's  theorem  may  be  applied  to 
show  that  the  requirement  is  satisfied  for  any  distribution  of  errors 
with  a  standard  deviation,  o-,  such  that 

k  (T  £  B 

C  =  P{|X-Xj  £k(r>>l  -  -V 

k"^ 

The  last  relation  may  be  solved  to  give  k  <  1/  tj  1  -  C ,  so  that  both 

equations  are  satisfied  if  we  have  er  £  B  n/  1  -  C.  More  precise  bounds 
on  cr  for  such  a  requirement  can  be  obtained  if  the  distribution  of  errors 
is  known.  For  example,  with  a  normal  distribution,  the  following 
table  applies  for  this  requirement; 

C  0. 5  0. 6  0. 7  ,  0. 8  0.  9 

cr£  1.481B  •1.187B  0.  965B  0.780B  0.  608B 


Variance  relations  for  the  errors  in  a  computed  quantity  in  terms 
of  the  variances  of  the  observed  quantities  have  been  discussed.  When 
the  variance  of  the  computed  quantity  is  prescribed,  limits  are  placed 
on  the  variances  of  the  observed  quantities.  The  allowable  range  of 
instrument  errors  could  be  apportioned  among  the  different  observed 
quantities  in  many  ways.  An  arbitrary  rule  that  has  been  proposed  to' 
resolve  this  indeterminate  problem  is  to  apportion  the  errors  such  that 
the  permitted  errors  in  the  observed  quantities  have  equal  effects  on 
the  errors  of  the  computed  quantities.  Thus,  if  X  is  the  computed 
quantity,  o^,  a^,  are  observed  quantities ,  and  the  variance 

equation  is 


2  a  2^.2,  ,  .  2 

cr„  =  A,  (T  +  A,  cr  +  .  .  .  +  A-.  tr 
X  1  Ujj  2  N 
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2  2  2  X 

choose  A  ,  tr  =  A-cr  =  .  .  .  =  A,,  cr  S  .  As  the  quavitities 

1  2  N  N 

A.  are  not  constant  but  are  functions  of  the  geometry,  a  conservative 
2  2 

rule  would  put  cr  ^o-^/NA.  .  The  extension  to  the  case  of  several 

u .  X  1  , 

1  Max 

dependent  variables  is  obvious. 

We  will  now  consider  weig?-.ting  factors  for  computed  quantities 
when  redundant  data  have  been  obtained.  We  shall  consider  the  case  in 
which  the  quantity  X  may  be  uniquely  specified  by  any  three  of  four 
measurable  quantities.  This  means  that  four  functional  relation¬ 
ships  exist  between  the  quantity  X  and  the  measured  quantities,  , 

Z^i  Z^)  and  Z_^.  We  shall  indicate  these  relationships  as 

X,  =  Xj  (Zj.  Zj.  Z^l 

Xj.X^tZj.  Z3,  Z^) 

X3-X3(Z,,  Z3.  Z^l 

X,  =  X^(Zj,  Z3,  Z3I 


If  small  errors  are  assumed, 


dX. 

1 


4 

y 

j 


9X. 

1 

JZT 

J 


dZ.- 

J 


I 

the  differential  errors  ia  X^  are  given  |by 


and  if  the  two  measured  quantities  Z^  and  Z^  are  uncorrelated  (see 
Supplement),  the  expression 


4  4 


m  T^i  i 


8X.  8X.  _ 

8^  ^TT 

m  i 


may  be  reduced  to 


_  ox. 

'"''j  "  z  ^ 


i 


8X.  8X. 
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The  computed  quantities  X^(i  =1,  2,  3,  4)  are  now  to  be 
weighted  to  get  the  best  estimate  of  X.  Thus,  if 
4 


i=  i 


where  is  the  value  of  the  weight  for  the  computed  X^,  then  the 
quantity 

_  ‘  4  4 

(dX)^  =  y  y  W.W.  dX.  dX. 

21  A/  1  J  1  J 

i=l  j=l 

must  be  minimized  with  respect  to  the  weighting  factors,  subject  to  the 
constraint 

i=  1 

If  the  undetermined  multiplier,  is  introduced,  this  constraint  may 
be  written  as 


2X 


4 

(l  -  .  0 

i=  1 


so  that 


4  4 


(dX)"  =  Z  Z  (‘  •  Z  '"i) 

i=lj=l  i=l 


2 

To  minimize  (dX)  with  respect  to  the  weighting  factors,  W^,  we  equate 


a(dX)'^ 

aw. 

1 


,  (i  =  1 ,  2  ,  3,  4),  to  zero.  Thus 
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a  (dx) 
aw. 


2 


4 


W.  dX.  dX.  +  N.  =  0, 
J  1  J 


(i  =  1.  2,  3.  4) 


These  four  equations  in  the  weighting  factors  have  a  solution 


W.  =  XB.(dXj^.  dX^.  dX3.  dX^),  (i  =  1,  2,  3,  4) 

where  the  are  obtained  from  diagonalizing  the  dX^^  dXj  matrix. 

Returning  to  the  normalization  constraint 

^  W.=\(B^  +B2  +  B3+B^)  =  1 
i  =  1 

the  value  of  \  is  found  to  be 


i  =  1 

and  finally  the  weighting  factors  are  given  by 
B. 

W.  =  (i  =  1.  2.  3,  4} 

j=  i 

Applications  and  Conclusions 


As  an  application  of  the  equations  derived  in  the  geometry 
portion  of  this  report  and  as  an  application  of  the  method  for  weighting 
the  results,  let  us  consider  the  following  example.  Assume  a  two- 
station  situation  in  which  the  base  line  f  is  known  exactly;  the  azimuth, 
elevation,  and  range  of  the  point  P  are  measured  from  one  station  (call 
these  Aj^,  E^,  and  pj^)  while  only  the  azimuth  and  elevation  of  the 

point  are  measured  from  the  second  station  (call  these  A2  and  E^).  It 

is  then  possible  to  apply  the  two-station  formulas  given  in  Cases  I 
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and  III  above.  Only  the  expressions  relating  range  to  the  measured 
quantities  will  be  used  in  this  example.  We  have  from  Case  I, 

cos^E^  sin^(A2  -  A^  )  +  4  sin^A^  -  cos  Aj^  sin  A^cos^E^^  sin(A2-Aj^)| 

R  J  - 2 - 2 - - 

4  cos  Ej  sin^(A2  -  A^^)  I 


and  from  Case  III 


'^2=|cJ 


^(-1 


^  +  cos^A^cos^E^^ 


+f  cos  A^cosE^ 


M'- 


i  ( 1  -  cos  A^cos 


From  these  two  expressions  it  is  desirable  to  find  weighting 
factors  which  give  a  best  estimate  of  R,  assuming  that  the  errors  (or 
more  exactly,  the  variances)  for  the  measured  quantities  are  known. 
That  is,  we  wish  to  find 


R  =  W^R^  + 

such  that  (dR)^  is  a  minimum,  and  subject  to  the  constraint 
Wi  +  W^  =  1 

In  a  straightforward  manner,  the  differential  errors  in  R^  and  R^ 
are  found  to  be 

dR^  =  dA^  +  dA^  +  dE^ 

and 

dR2  =  Lj  'iPj  +  Lg  ’^■^2  ^  ^3  *^^2 


where 


= 


KdA^J 


M 


sinA^ 


2  sin  cos(A2  -  A^) 


4  sin(A2  -  Aj^)  cos^E^  sin^(A2  -  A^) 


sin  A^  cos  Aj  cos  (A^  -  Aj^) 
2 

4  sin  (A^  -  Aj) 
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K. 


=  A)  = 

\8A,  j 


M 


2  sin  A^  cos  A^ 


cos  A^  cos  Aj 


2  2 

cos  Ej^  sin  (A^-Aj^)  sin(A2-Aj^)4 


.  2 


2  sin  A^  cos(A2  -  A )  sin  A^  cos  A^^  cos(A2  -  A ^ ) 

— z - ; — 3 -  - z - 

cos  Ej  sin  (A^-Aj)  4  sin  (A^-Aj^) 


K, 


X 

VSEj,/'  " 


M 


2 

2  sin  A^  sin  E 

2  3 

sin  (A^  -  Aj^)  cos  E^ 


With 


M 


•  2. 

,  sin  A- 

_  + _ £_ 

4^  .._2^  .,.2, 


sin  A^  cos  Aj^ 


cos  E.,  sin  (A,  -  A . ) 


4  sin(A2  -  A^) 


and 


- 

"  \  8p^  ! 


N 


2p,  + 


i  p  ^  cos  A^  cos  E^ 


^  Fz  ,  2  ,  ,  2~  2„  \ 

Ip,  -i  ’  1  -  cos  A-,  cos  E- I 

i  2  2/ 


/  ®^2  \ 

^2=  =  ^ 


sinA^  cos  E- 


-cos‘'A2Cos‘'E2  [  ^  1 


cos  A^  cos  E^ 


fz  TzTi  F  2“ 

/pj  -i  (1  -  cos  A^  cos  E_ 
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=  A)" 

V8E,  ) 


N 


-i  sin  cos  J  cos^A^cos^E^^l^l 


I  cos  A^  cos  E^ 


^  2 


JpJ-i^{l-cosA2Cos2E^ 


with 


N 


=  H 


2  ,  «2  /  3 

+ i  ( - ^ 


Z  2  \ 

,  ^  +  cos  A,  cos  E^  1 

\  4  2  Z) 


+  I  cos  A^  cos  E^ 


2t:.  \ 

COS  cos  E-  } 
Z  Zj 


r 


From  the  previous  development,  we  take  the  derivative  of 


=  1 1  *i'*j  (‘  -  E  ) 


1  J 


=  Wj  (dR^r  +  2Wj^W2  dR^  dR^  +  (dR^)  -  2\(1  -  W.^  -W^) 

=  wf  (kA^  +  ^  +  2W,W,(k,L,  0-^  ^ 

1  V  1  A^  2  A^  3  Ej/  1  2V  2  2  A^/ 

+  +L^o-^  +  )  -  2\  (l  -W,  -  W,^ 

2\lp^  2  ^2  ^  ^2^  \  1  2/ 

with  respect  to  the  weighting  factors,  to  obtain  the  two  equations 
I'Z  Z.  ^  1,2  2  j^2  2  '\  A 

■Krr  4-KrT-  +K«-r  >W  4-*F 


/l^T  ^  ■'iW  XT^^  XT'- 

J,K^L-,cr.  tW,  +  *;L,(r  +  L^t*  +  L^c 

V,  2-  2  A^/  1  \  lp.^  2A 


K^L,o-^  ^ 

W,  = 

-  X. 

2  2  A^,' 

2 

+  L^o-^  ^ 

= 

3  E^; 

2 
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Through  use  of  the  normalization  condition,  these  equations  may  be 
solved  to  give 


W 


1 


W 


2 


D 

J_ 

D 


♦  (‘-2  -  K2^2)4,  + 

"  (k2  - 


where 


D 


- 


\  2  2 
'^2-^2)  <^A 


( 


^  „2  2  ,,22 
+  K,cr_  +  L  ,  0- 

2  2  ^1  ^  P 


-u  T  2  2 

+  L,o-„ 

1  P  ^2 


To  carry  the  example  further  the  calculation  of  numerical  values 
for  the  variances  (or  average  squares  of  errors)  of  Rj^,  R^,  and  the 

weighted  value  of  R  will  be  outlined.  The  values  of  the  parameters 
chosen  have  no  special  significance.  Values  for  I,  E^,  and  Aj 

were  assumed,  and  values  of  E^  and  A^  were  then  computed  to  give 

"near"  closure  of  the  various  triangles.  These  values  are 

I  =  100  feet 

=  2000  feet 

Ej  =  20° 

A^  =  60° 

E^  =  20°29' 

A^  =  62°42' 

For  this  particular  choice  of  parameters,  it  is  found  that 
Rj  =  1984.  69  feet 
R^  =  1980.79  feet 

These  values  are  in  very  good  agreement.  If  the  values  of  and  A^ 

had  not  been  chosen  to  give  "near"  closure  of  the  triangles,  a  much 
larger  disagreement  between  Rj  and  R^  would  have  resulted. 
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The  other  computed  quantities  of  interest  from  this  set  of 
parameters  are 

=  4.  2610  X  10^  feet 

=  -4.  1438  X  10“^  feet 

=  7.  3916  X  10^  feet 

=  9. 9884  X  10"^ 

=  4.  0197  X  10^  feet 

=  1-  5828  X  10^  feet 

The  variances  for  the  distance  and  angular  measurements  are  then 
assumed  to  be 

<r^(pj^)  =  lO^ft^ 

tr^  (all  angles)  =  9X  10”^ 

This  gives  the  variances  for  Rj^  and  as 

0-0  =  31799.  ft^ 

^1 

0-0  =  99.7  86  ft^ 

2 

By  the  process  outlined  above,  the  weighting  parameters  are  found  to 
have  the  values 

=  .00359 

=  .99641 

The  weighted  value  of  R, 

R  =  WjRj  +  W^R^ 

is  then  found  to  be 

R  =  1980.  8  ft 


with  a  variance  of 
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0--^=  99.  3733 

By  direct  comparison  then,  it  is  seen  that 

2  ^  2  ^  2 

<  ‘"R2<  '^Rj 

which  points  out  the  value  of  the  weighted  average  for  the  computed 
quantity,  in  this  case  R. 

For  emphasis,  it  is  pointed  out  again  that  when  a  particular 
measuring  scheme  has  been  set  up  and  the  variances  have  been 
established  for  the  measured  quantities,  then  the  variances  of  the 
computed  quantities  may  be  calculated.  Or,  if  the  maximum  allowable 
variance  is  specified  for  a  particular  computed  quantity,  limits  may  be 
established  for  the  measured  quantities;  i.  e.  ,  the  maximum  errors 
that  the  measuring  devices  could  have  for  the  error  of  the  computed 
quantity  to  be  less  than  some  given  amount  could  be  established.  To 
determine  those  limits  for  the  errors  of  the  measuring  devices,  however, 
the  functional  relations  given  previously  between  the  computed  and 
measured  quantities  and  the  differentials  of  these  functional  relations 
are  required.  The  method  of  weighting  parameters  is  required  only 
when  a  redundancy  of  data  occurs. 
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SUPPLEMENT 

PROPAGATION  OF  COVARIANCE  OF  ERRORS 


Suppose  certain  variables  y-j^,  k  =  1,  ?,,  .  .  .  ,  m,  functions  of 
another  set  of  variables  x^,  s  =  1,  2,  .  .  .  ,  r,  are  given  as  functional 
relations  =  yj^  (x^,  x^,  .  .  .  ,  ^j.)-  Linear  relations  for  small  errors 
in  y  due  to  small  errors  in  x  are  then  given  by  the  differentials, 


(1) 


These  errors  may  be  viewed  in  a  statistical  sense  if  the  linear 
error  relations  (1)  hold  for  the  great  bulk  of  the  errors,  so  that  these 
relations  may  be  averaged  over  the  errors.  If  all  the  systematic 

errors  and  biases  have  been  removed  so  that  the  errors  dy,  and  dx. 

^  J 

appearing  in  (1)  are  accidental  errors,  it  may  be  inferred  that  the 
average  errors  are  zero  for  both  x  and  y. 

The  addition  of  a  second  subscript  t  to  the  quantities  occurring  in 
equation  (1)  may  be  made  to  indicate  the  results  of  the  t-th  experiment 
in  measuring  the  errors. 


II 

>> 

V  1 

1 

(2) 

'^yst  = 

% 

1 

(3) 

j 


The  covariance  of  the  errors  in  y  may  be  obtained  by  taking  the 
product  of  equations  (2)  and  (3),  and  averaging  over  all  the  experiments. 


'^^kt  '^^st 


I. 


dx...  dx... 
It  jt 


(4) 


Lim 
h-»  00 


h 


h 

) 

t=l 


^^kt  '^^st 


Lim  1 
h—  00 


h 

\'  Y 
Z.-  Zv 
t=l  i,j 


^^k  ^^s 

Bx.  8x. 
1  J 


dx..  dx.. 
it  jt 


(5) 
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or 


or 


i-j 

i.  j 


9x. 

9x.  '^x.x. 

1 

J  1  J 

(6) 


Equation 
errors .  If  the 


(6)  is  a  relation  of  the  propagation  of  covariance  of 
errors  dx^  and  dx.  are  independent  for  i  j,  then 


cr  reduces  to  6'?,  where  S'?  is  the  Kronecker  delta,  and 
x.x.  X.  1  1 

1  J  1 

equation  (6)  reduces  to 


9yc 


dx.  } 
1  ' 


(7) 


This  gives  the  covariance  of  errors  in  y  in  the  event  the  errors  in  x 

are  independent.  In  particular,  for  y.  =  y  .  there  results  for  this 

s 

case  of  independent  errors  in  x. 


L  /' 


1 


(8) 


Equations  (6),  (7),  and  (8)  may  be  useful  for  interpreting  in  a 
statistical  sense  the  errors  in  functions  of  a  set  of  variables,  which 
may  be  either  dependent  or  independent. 
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APPENDIX  5 

LEAST-SQUARES  ADJUSTMENT  OF  TWO-STATION  ANGLE-ONLY 
POSITION  FIXES  WITH  RELIABILITIES  OF  ADJUSTED  DATA 


Introduction 


One  phase  of  the  present  program  of  investigating  scoring 
concepts  for  weapon -target  encounters  is  concerned  with  the  limitations 
imposed  by  errors  in  the  instrumentation.  Computational  methods 
making  use  of  redundant  data  can  improve  the  reliability  and  usefulness, 
of  the  data  taken,  and  this  paper  is  concerned  with  the  least-squares 
method  of  adjusting  redundant  data  and  with  the  reliability  of  adjusted 
data  and  of  quantities  computed  from  adjusted  data. 

Certain  of  the  scoring  concepts  being  investigated  are  based  on 
angle -only  observations  made  from  two  or  more  stations  separated  by 
base  lines  of  known  lengths  and  orientations.  At  each  such  station 
observations  are  made  of  the  azimuth  and  elevation  angles  of  the  line  of 
sight  to  an  object  whose  position  is  to  be  determined.  It  is  clear  that 
if  simultaneous  azimuth  and  elevation  data  are  available  from  each  of 
two  or  more  stations,  the  data  are  more  than  the  minimum  needed  to 
determine  the  observed  object's  position,  and  there  is,  therefore, 
redundancy  in  the  data. 

The  method  of  least  squares  applied  to  redundant  angular  data 
gives  adjusted  observations  with  minimum  variances.  It  was  desired 
to  know  more  about  the  extent  of  the  improvement  in  reliability  obtained 
by  adjusting  the  data  and  about  the  least-squares  process  itself  as 
applied  to  this  particular  problem.  Since  a  survey  of  the  literature 
did  not  obtain  any  useful  information,  a  limited  investigation  was  under¬ 
taken.  This  paper  is  concerned  with  setting  up  a  least-squares  process 
to  adjust  the  azimuth  and  elevation  angles  of  an  object  observed 
simultaneously  from  two  stations  separated  by  a  base  line  of  known 
length.  The  angular  data  are  assumed  to  have  equal  variances  and  to 
be  uncorrelated  (independent)  and  unbiased.  Expressions  are  obtained 
for  the  separation  of  the  object  from  either  of  these  observation 
stations.  Expressions  for  the  variances  of  the  adjusted  data  and  of  the 
computed  ranges  are  also  developed.  It  is  shown  that  least -square 
adjustment  of  the  data  reduces  the  probable  errors  in  both  measured 
and  computed  quantities  by  about  50  percent  on  the  average. 

One  advantage  of  this  process  over  certain  other  least-squares 
procedures  is  that  preliminary  estimates  of  position  are  not  required. 
The  procedures  used  to  adjust  the  data  and  to  obtain  variance  estimates 
for  the  adjusted  observations  and  computed  quantities  are  straight¬ 
forward  and  fairly  simple.  A  disadvantage  is  that  the  one  condition 
equation  is  required  to  force  the  adjusted  rays  from  two  stations  to 
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intersect,  and  an  additional  condition  equation  would  have  to  be  used 
for  each  added  station;  as  a  result,  the  adjusted  data  have  non-zero 
covariances.  The  process  could  be  readily  generalized  to  include 
correlated  data  of  differing  reliabilities  and  to  include  more  stations 
and  more  observational  data;  the  process  is  not  limited  to  angular 
data.  For  the  present  purpose  of  determining  the  variances  of  the 
adjusted  data  and  position  fixes,  however,  these  generalizations  were 
not  needed  and  were  not  developed. 

A  Condition  Equation,  Necessary  and  Sufficient  for  Rays  from  Two 
Stations  to  Intersect 


The  elevation  and  azimuth  angular  data  are  assumed  to  have  been 
corrected  for  all  known  errors.  These  include  biases,  instrument 
errors,  refraction  errors,  etc.  Since  some  random  error  will 
normally  remain  uncorrected,  the  two  rays  (lines  of  sight)  defined  by 
the  corrected  angular  data  will  seldom  intersect.  The  data  are  to  be 
adjusted  by  the  method  of  least  squares  to  give  unbiased  mean  values  of 
minimum  variances,  and  an  additional  restricting  condition  is  made 
that  the  adjusted  rays  intersect. 

Figure  1  is  a  sketch  of  the  geometry  of  the  problem.  A  Cartesian 
coordinate  system  is  shown;  the  distance  between  the  observation 
stations  at  O  and  Q  on  the  Y  axis  is  denoted  by  I  .  The  adjusted  position 
of  the  observed  point  is  at  P  where  the  adjusted  rays  intersect. 


Figure  1 

Adjusted  Two-Station  Geometry 
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From  the  geometry, 

Z  =  Sj^  tan  Ej^  =  tan  E2  (1) 

By  the  law  of  sines, 

i  sin  I  sin  A. 

S  =  - i-,  S  = - i-  (2) 

sin  (A^  -  Aj^)  sin  (A^  -  Aj) 

Substitution  of  equations  (2)  into  equation  (1)  gives  the  condition 
equation  for  the  adjusted  observations, 

f  =  sin  A^  tan  E^  -  sin  A^  tan  E^  =  0  (3) 

Equation  (3)  may  also  be  obtained  by  equating  the  X/Z  ratios  of 
the  adjusted  observed  point  as  measured  at  the  two  stations.  The 
foregoing  development  assumes  that  the  algebraic  signs  of  A^,  A^, 

and  (A^  -  Aj^)  are  the  same;  E^  and  E^  should  also  have  the  same  signs. 

Equation  (3)  may  be  shown  to  be  a  necessary  and  sufficient 
condition  for  the  adjusted  rays  to  intersect,  possibly  at  infinity.  The 
necessity  is  obvious  from  the  foregoing  development.  If  more  than 
two  stations  were  used,  an  additional  condition  equation  would  be 
needed  for  each  such  additional  station.  Since  these  additional  restric¬ 
tions  are  not  needed  for  the  present  purposes,  they  will  not  be 
developed. 

The  Least-Squares  Adjustment 


Suppose  simultaneous  azimuth  and  elevation  angles  A^^,  Ej^,  A^, 

and  E^  are  measured  at  two  observation  stations  as  shown  in  Figure  1, 

In  general,  equation  (3)  is  not  satisfied  exactly,  so  that  a  small  error 
remains, 

f(Aj,  E^,  A^.  E^)  =  d  (4) 


Using  the  linear  terms  of  a  Taylor's  series  expansion  of  f,  and 
representing  the  changes  to  be  added  to  the  raw  data  to  give  the 
adjusted  data  as  6Aj^,  6Ej^,  etc.,  the  linearized  condition  equation  for 

the  adjusted  observations  is  obtained: 


f  =  d  + 


9f 

8A, 


6A,  + 


ar 

9E, 


6E,  + 


9f 

9A, 


6  A. 


9f 

9E_ 


6E2  =  0 


(5) 
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Upon  evaluating  these  partial  derivatives  at  the  known  approximate 
values  of  the  angles,  equation  (5)  becomes 

f  =  d  +  a^  6Aj^  +  a^  6Ej  +  a^  6A2  +  a^  6E2  =  0  (6) 

in  which 

=  -  cos  A^  tan  E2 

■  A  2t:' 

a  =  sin  A  sec  E, 

221  (yj 

a^  =  cos  A2  tan  Ej 
2 

a,  =  -  sin  A,  sec  E^ 

4  12 


The  adjusted  observations  are  related  to  the  observed  data  as  follows: 


A,  .  ,  =  A,  +  6A, ,  etc. 

1  adjusted  1  1 


(8) 


In  most  cases  of  practical  interest,  the  errors  in  A^,  Ej,  A2»  E2, 

are  uncorrelated  and  may  be  considered  to  have  equal  variances,  (Tq. 

In  such  cases,  the  observations  all  have  equal  weight,  and  the  least- 
square  problem  is  to  minimize 

e^A^  +  e^Ej  +62a2  +  62e2, 

subject  to  the  linearized  condition  equation  (6). 

The  easiest  way  to  proceed  is  to  make  use  of  a  Lagrangian 
multiplier,  X.,  and  to  minimize 

H  =  6^Aj  +  6^E^  +  +  6^E2  -  2kf  (9) 

Upon  taking  the  partial  derivatives  of  H  with  respect  to  A^  Ep  etc.  , 
and  equating  them  to  zero,  the  following  normal  equations  result; 
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26A^  -  2Ka^  =  0 

26E  -  2Xa  =  0 

(10) 

26A2  -  2\a2  =  0 

26E.  -  2Xa.  =  0 

2  4 

The  normal  equations  (10),  together  with  equation  (6),  may  be 
solved  for  the  unknowns,,  with  the  result 

6A^  =  -aj^d/(a^  +  a^  +  a^  +  a^) 

2  2  2  2 
6E,  =  -a^d/(a.  +  a^  +  a,  +  a.) 
i  ^  1  ^  J  4 

6A2  =  -a2d/(a^  +  a^  +  a^  +  a^) 

6E2  =  -a4d/(a^  +  +  a^  +  a^) 


These  are  the  corrections  to  be  applied  to  the  observed  data 
according  to  equation  (8).  They  should  satisfy  the  linearized  condition 
equation  (6),  which  serves  as  a  check  for  small  corrections.  If  the 
corrections  should  be  too  large  for  equation  (6),  obtained  from  linear 
terms  in  Taylor's  series,  to  be  applicable,  an  iteration  process  could 
be  used  to  obtain  a  better  answer. 

Reliability  of  Adjusted  Data  and  Computed  Functions 

Considerable  interest  may  be  attached  to  the  reliabilities  of  the 
adjusted  data  and  those  of  functions  computed  from  these  data.  In  the 

2 

following  it  will  be  presumed,  that  the  variance,  cTq,  of  the  unadjusted 

angular  measurements  is  known  in  advance.  However,  it  will  be  seen 
that  this  value  may  be  estimated  from  the  data. 

The  adjustments  to  be  made  to  the  data  satisfy  equation  (6).  This 
equation  may  be  written  in  the  form 


-d  =  aj^6Aj  +  a^dEj^  +  a.^6A^  +  a^SE^ 


(12) 


Upon  squaring  this  equation  and  averaging,  there  results 


Average 


,  2  ^  2  ^  2  ^  2.  2 

=  (a^  +  a^  +  a3  + 


(13) 
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This  result  may  be  used  with  equation  (11)  to  obtain  the  variances 
of  the  adjusted  observations.  Upon  squaring  equations  (11)  and  averag¬ 
ing,  making  use  of  equation  (13),  one  obtains 


) 


Adjusted 


2  2  ,,  2  2  2  2 , 
ai  ^Q/(a^  +  a2i-a3+a4) 


,2.  22,,2,2^2,^2. 

1  Adjusted 


,2,  22.,2,2,22v 

*^^7  A.I-  I-  A  ~  ^3  ^2  ^3  ■'■^4^ 

2  Adjusted 


(14) 


I  2  > 

2  Adjusted 


2  2  ,  ,  2  2  2  2. 

”  ^4  °’0'^  ^"2  '^3  ^4^ 


These  results  show  that  on  the  average  the  variance  of  any  one  of 
these  adjusted  observations  is  only  one  fourth  of  the  variance  of  the 
unadjusted  observation,  since 


+  0-; 


+  tr , 


2^  Adjusted 


2 

•"0 


while 


.U  2  ,  2  . 

+  0"a  ' 

2  2  Unadjusted 


(15) 


It  may  be  noted  from  equation  (14)  that  as  any  coefficient  a.^  tends 
to  zero  (while  at  least  one  of  the  other  0,  j  i),  the  variance  of 
the  corresponding  adjusted  observed  variable  tends  to  zero. 


A  result  of  placing  the  restricting  condition,  equation  (6),  on  the 
adjusted  variables  is  that  they  no  longer  are  independent,  but  have 
covariances  which  are  in  general  different  from  zero.  These  covari¬ 
ances  may  be  readily  found,  using  equations  (11)  and  (13).  Upon 
taking  the  products  of  pairs  of  equations  listed  under  (11),  averaging, 
making  use  of  equation  (13),  there  results 
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'A„E 


11 


^1^2 


^1^2 


'E,  A 


r‘-2 


’E,E 


12 


■^2^2 


^1^2'^0^  ^^1  ■*■  ^2  ■*■  ^3  ■'■  ^4^ 

^I^S'^q/^^I  +  ^2+  ^3+ 

2.,  2^  2^  2^  2. 

^1  ^2  ^  ^3  ^  ^4 

a.^a,^cr^/ +  a^,  +  a^  +  a^) 

^2^4‘"o/(^l  +  ^2  +  ^3  + 

2  ,  ,  2  2  2  2, 

^3^4°'0^  ^^1  ■*■  ^2  +  ^3  ■*■  ^4) 


(16) 


This  dependence  of  the  adjusted  observations  is  sometimes  over¬ 
looked  in  variance  computations.  The  variance  of  any  function  of  the 
variables  may  be  computed  from  their  adjusted  variances,  but  the 
covariance  terms  should  be  included  in  such  a  computation.  The  covar¬ 
iance  terms  are  usually  smaller  than  the  variance  terrns  and  there  may 
be  cancellation  effects  due  to  algebraic  signs.  One  has  always 
cr„  Y  =  o-vCr.,,  p,  with  a  correlation  coefficient  p  such  that  (pi  <1. 

^1^2  ^1^2 

2 

It  is  interesting  to  note  that  the  same  average  value  of  d  is 
obtained  by  using  either  the  unadjusted  observations  or  the  adjusted 
observations;  this  serves  as  a  check  on  the  calculations.  It  may  be 
noted  that  the  coefficients  a^^  are  simultaneously  equal  to  zero  only  if 

the  observed  point  li^s  along  the  Y  axis.  In  this  case  no  change  in  the 
observations  is  made  by  the  least-squares  adjustment.  Also,  the  two 
rays  intersect  at  an  indeterminate  point  along  the  Y  axis,  and,  as  will 
be  seen  from  later  equations,  the  variance  of  the  range  is  infinite. 

This  indicates  that  one  should  strive  to  have  the  observed  point  not 
located  along  the  Y  axis  or  even  close  to  it.  ; 

Expres  sions  ,for  the  range  may  be  obtained  by  the  use  of  trig¬ 
onometry.  By  use  of  the  law  of  sines  for  triangle  POQ  in  Figure  1, 
one  obtains 
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= - 

cos  A2COS 


’-I- 


£j  1  -  cos^E^  cos^A2 


(17) 


I  2  2  I  2  2 

,  /I  -  cos  Aj^cos  E^+  cosA.jCosEj^/1-cos  A^cos  I 
Alternate  forms  are  less  complies 


with  a  similar  expression  for  R^. 

=  i  sin  A^/cos  Ej^  sin  {A^  -  Aj^) 

A,) 


R 


i  sin  A  j/cos  E^ 


E^  sin  (A^  -  jT-j^ 


ing 


2 

equations  (18)  may  be  taken,  with  the  foil 


The  differentials  of 
results; 

dR^  =  Mj,dA^  +  ^3'^-^2 

dR^  =  Nj^dAj^  +  N^dE^  +  N^dA^ 


where 


Mj 

=  i 

sin 

CSC  (A^ 

-^1 

cot  (A2 

-  A3)  sec 

^1 

M2 

=  t 

sin 

^2 

CSC  (A^ 

-  A3) 

sec  E3 

tan  E3 

M3 

=  t 

sec 

^1 

CSC  (A^ 

-Ai) 

[  COS  A 

2  -  "^"^2 

cot 

(A2 

■  A3)] 

^1 

=  i 

sec 

^2 

CSC  (A^ 

-  Ai) 

[  COS  A 

3  +  sin  A3 

cot 

{A^ 

-A3)] 

^2 

=  i 

sin 

Ai 

CSC  (A^ 

-A3) 

sec  E2 

tan  E2 

N3 

= 

f  sin  A 

sec  E2 

CSC  (A2  -  A3 

)  cot  (A2  ■ 

) 

result: 


Upon  squaring  and  averaging  equations  19,  the  following  variances 


(4 ) 


1  Adjusted 


(21) 


U  "'•3“A^A2 


2  3' 


^EiA2  J 


Adjusted 
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(4 ) 

2  Adjusted 


^,2  2  »,2  2  ,  »,2  2 

^I^'Aj  ■*■  ’^2°'E2  ■*■  ^^1^2‘^A^E, 


+  .  +  ^N^NjO-g  .  J 

1  2  2  2J  Adjusted 


Since  the  range  is  computed  from  the  adjusted  variables ,  its 
variance  will  be  approximately  proportional  to  the  variances  of  the 
adjusted  variables.  The  variance  of  the  range  is  thus  on  the  average 
approximately  one  fourth  of  what  it  would  be  if  the  range  were  computed 
from  the  unadjusted  variables;  the  standard  error  in  range  is  reduced 
by  a  factor  of  about  one  half.  Cartesian  coordinates  of  the  observed 
object  could  of  course  be  computed  from  the  snherical-polar  coordi¬ 
nates,  and  their  variances  could  be  readily  obtained  by  the  use  of 
differentials . 
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APPENDIX  6 

APPLICATION  OF  POLYNOMIAL -BASED  SMOOTHING  AND 
INTERPOLATING  FORMULAS  TO  SCORING  PROBLEMS 

Many  formulas  are  available  for  smoothing  and  interpolating  data 
taken  at  a  constant  time  interval.  Most  of  these  formulas  are  based 
on  the  premise  that  the  function  of  the  time  represented  by  the  data  on 
which  these  operations  are  carried  out  is  given  with  sufficient  accuracy 
over  a  short  interval  by  a  polynomial  of  a  specified  low  degree.  Since 
the  results  will  not  be  satisfactory  in  general  unless  this  assumption 
is  valid,  it  is  desirable  to  determine  the  true  nature  of  the  function 
represented  by  the  data.  The  purpose  of  this  paper  is  to  show  that 
polynomial-based  formulas  are  suitable  for  use  on  certain  types  of 
data  obtained  for  scoring  the  attack  of  a  missile  on  a  target. 

With  trajectories  of  long  duration,  it  may  be  difficult  to  specify 
a  priori  the  probable  character  of  a  variable  function  of  the  time. 
However,  with  trajectories  of  short  duration,  where  accelerations 
have  limited  time  to  affect  the  variables,  the  probable  behavior  of  a 
variable  may  often  be  inferred.  Thus  a  linear  relative  trajectory 
provides  a  good  approximation  for  a  relative  trajectory  of  sufficiently 
short  duration.  By  a  linear  relative  trajectory  is  meant  a  trajectory 
of  the  weapon  relative  to  the  target  which  is  given  by  a  linear  vector 
function  of  the  time, 

R  =  R  +  R  t  (1) 

o  o  '  ' 

* 

with  R  and  R  constant  vectors.  The  addition  of  a  small  non-zero 

o  o 

acceleration  term  to  the  right  side  of  equation  (1)  would  not  affect  the 
conclusion  to  be  drawn  in  the  following  discussion. 

The  components  of  equation  (1)  in  a  rectangular  Cartesian 
coordinate  system  are  as  follows: 

X  =  X  +  X  t 
o  o 

Y  =  Y  +  Y  t  (2) 

o  o 

Z  =  Z  +  Z  t 
o  o 

For  short-duration  trajectories  there  is  considerable  probability  that 
the  relative  Cartesian  coordinates  of  the  weapon  with  respect  to  the 
target  are  well  approximated  by  low-degree  polynomials  in  the  time. 
Polynomial -based  smoothing  and  interpolating  formulas  should  thus  be 
suitable  for  application  to  the  relative  Cartesian  coordinates  X,Y,Z. 
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The  case  is  different  for  relative  polar  coordinates,  as  will  be 
seen.  Let  the  polar  coordinates ,  ra,nge,  azimuth,  and  elevation,  be 
considered  as  functions  of  the  time,  considering  first  the  range,  R. 

=  R  •  R  s  R^  +  2R  •  R  t  +  R^t^  (3) 

O  O  O  0 

7  ~  7  - 

|r|  =  {r'^)2  s  (a  +  bt  +  ct  (4) 

Equations  (3)  and  (4)  show  that  R^  for  a  linear  relative  trajectory 
may  be  well  approximated  by  a  polynomial  of  low  degree  while  R  may 
not,  in  general.  This  suggests  that  better  results  should  be  obtained 
2 

if  R  raither  than  R  is  smoothed  or  interpolated.  If  R  is  smoothed  or 

interpolated,  good  results  should  be  obtained  only  for  R  and  R  closely 
aligned.  °  ° 

The  figure  shows  the  non-rotating  Cartesian  coordinate  system 
X,Y,Z  with  origin  O  at  the  target.  The  weapon  is  at  P,  and  it  will  be 
seen  that  the  azimuth  and  elevation  angles  are  given  by 

d  +  et 
f  +  gt 

(5) 

h  +  kt 

(a  +  bt  +  ct^)^ 


.  A  X 

tan  A  =  ^  ~ 
sin  E  =  ~ 


X 
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Thus  A  and  E  are  not  given  even  approximately  by  low-degree  poly¬ 
nomials  in  t  over  any  appreciable  time  interval.  If  the  variables  A 
and  E  are  observed,  and  are  operated  on  with  polynomial -based 
smoothing  or  interpolation  formulas,  the  results  should  not  be  expected 
to  be  very  good  over  appreciable  time  intervals. 

The  angular  rates  A  and  E  of  R  may  be  determined  by  differentia- 

'  K  • 

tion  of  equations  (5).  It  is  found  that  A  - - y.  Thus,  1/A 

M  +  Nt  +  Qt"^ 

should  be  given  approximately  by  a  low  order  polynomial  in  t,  and  for 

1  • 

best  results  1/A  rather  than  A  should  be  used  with  polynomial-based 
smoothing  or  interpolating  formulas.  On  the  other  hand,  neither  E 
nor  1/E  appears  to  be  given  approximately  by  a  low-degree  polynomial 
in  t  over  any  appreciable  time  interval.  However,  for  R^  an"’  R^  some¬ 
what  closely  aligned,  the  quantity  (sec  E)/E  would  be  give  .  approxi¬ 
mately  by  a  low-degree  polynomial  in  t,  so  that  operating  on  this  quantity 
should  be  expected  to  yield  more  valid  results  than  operating  on  E  or 
•  2  • 

E  alone.  It  may  also  be  shown  that  (sec  E)/A  should  be  closely 
approximated  by  a  low-degree  polynomial  in  t  for  linear  relative  tra- 

4 

jectories,  even  with  R^  and  R^  not  closely  aligned. 

In  conclusion,  it  appears  probable  that  relative-Cartesian- 
coordinate  data  for  a  short- duration  weapon-target  encounter  may  be 
well  adapted  to  numerical  methods  of  extrapolating,  interpolating, 
and  smoothing,  whereas  relative-polar-coordinate  data  does  not 
appear  to  be  so  well  adapted.  A  reasonable  procedure  for  performing 
these  operations  on  the  relative  polar  coordinates  would  be  to  first 
convert  the  data  to  Cartesian  coordinates,  perform  the  operations 
required,  and  then  convert  the  resultant  values  to  polar  coordinates. 
Certain  functions  of  the  relative  polar  coordinates  have  been  shown  to 
be  well  approximated  by  polynomials  (in  the  time)  of  low  degree  and  are 
therefore  well  adapted  to  the  numerical  operations  of  smoothing, 
interpolation,  and  extrapolation. 
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